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.simulated  in  a  realist  i<  manner.  Also,  the  results  ol  simulations  have  been  ’ 
compared  to  measurements  from  the  ESMR,  SCAMS,  and  H I KS  sensors  on  the 
Nimbus  VI  satellite,  which  was  the  most  complete  set  ol  microwave  and  infrared 
daLa  available  at  the  Lime  ol  this  study.  The  first  area  is  concerned  with 
the  investigation  ol  the  effect  of  clouds  and  precipitation  on  the  Nimbus  VI 
ESMR  channel  and  Lhe  use  ol  KSMR  vertical  and  horizontal  polarization 
components  for  the  determination  ot  the  liquid  water  content  and  cloud 
thickness.  In  conjunction  with  these  objectives,  we  have  developed  a  micro- 
wave  radiative  transfer  program  taking  into  account  the  polarization  effect. 
Utilizing  the  microwave  radiative  transfer  program,  comparisons  between  the 
computed  and  observed  brightness  temperatures  using  the  available  radiosonde 
and  surface  observations  are  first  carried  out  and  theoretical  sensitivity 
analyses  then  follow.  A  simple  emp i r ica 1 -theore t i ca 1  approach  is  subsequently 
taken  to  retrieve  the  cloud  liquid  water  content  and  thickness . -^The  second 
area  deals  with  the  use  of  both  infrared  and  microwave  channels  for  the 
cloud  parameter  inference  parameter i zat i ons .  The  infrared  and  microwave 
radiative,  transfer  programs  are  successively  developed  for  retrieval  of  the 
high  and  low  cloud  top  heights  and  thicknesses  and  the  cloud  liquid  water 
content  from  selected  I1IRS  and  SCAMS  channels  of  the  Nimbus  VI  satellite. 

Random  error  analyses  revealed  that  the  cloud  sounding  technique  employing  a 
combi  rat i on  ol  infrared  and  microwave  frequencies  appears  to  be  theoretically 
rigorous  and  practically  workable.  Thus,  applications  of  the  cloud  sounding 
technique  to  the  Nimbus  VI  11 1 RS  and  SCAMS  data  are  carried  out  in  which  we 
have  selected  a  number  ot  cases  associated  with  summertime  convective  cloud 
svstems  and  wintertime  large  scale  synoptic  cyclones  described  in  a  previous 
report.  The  report  presented  detailed  discussions  on  the  surface  and  upper 
air  observations,  radar  echoes,  precipitation  and  other  relevant  synoptic 
features  which  were  utilized  to  cross-check  with  the  retrieved  cloud  parameters 
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SECTION  1 


INTRODUCTORY  REMARKS 

In  our  previous  research  work  concerning  the  development  of 
cloud  sounding  techniques  from  orbi ti ng  meteorological  sateMites 
(Feddes  and  Liou,  1978;  Liou  and  Duff,  1979),  it  was  pointed  out  that  a 
proper  combination  of  infrared  and  n.icrowave  measurements  could  result 
in  a  significant  improvement  for  the  vertical  and  horizontal  mapping  of 
the  cloud  compositions  and  structure.  In  an  effort  to  follow-up  this  idea 
and  to  provide  a  theoretically  rigorous  foundation  for  cloud  sounding  from 
satellites  based  on  the  theory  of  radiative  transfer  in  scattering  and 
absorbing  atmospheres,  we  have  proposed  a  program  to  integrate  radiative 
transfer,  remote  sensing,  and  synoptic  meteorology  to  furthe.  advance  the 
development  of  cloud  sounding  techniques.  The  specific  objectives  of  the 
present  research  project  under  the  sponsorship  of  the  Air  Force  Geophysics 
Laboratory  have  been  (a)  to  develop  a  radiation  model  including  polarization 
for  microwave  imagery  channels,  (b)  to  develop  techniques  for  parameterizi ng 
the  satellite  microwave  and  infrared  imagery  and  sounder  data  for  the 
inference  of  cloud  compositions  and  structure,  and  (c)  to  utilize  the 
available  ground  truth  data  with  varying  degrees  of  visibility,  thickness, 
and  coverage  for  both  cirri  form  and  water  clouds,  including  sources  from 
weather  station  reports,  radiosondes,  radar  summaries,  and  infrared  pictures 
to  cross  check  with  the  retrieved  cloud  parameters. 

Within  the  scope  of  these  objectives  and  subject  to  the  time 
constraint  and  funding  availability,  we  have  utilized  the  comprehensive 
infrared  and  microwave  data  set  from  the  Nimbus  VI  satellite  and  the 
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simultaneous  synoptic  reports,  some  of  which  are  available  through  our  previ¬ 
ous  research  projects  sponsored  by  the  Air  Force  Geophysics  Laboratory.  In 
this  report,  we  present  two  areas  of  research  both  of  which  are  related  to 
the  determination  of  cloud  parameters  from  meteorological  satellites. 

In  Section  2  of  this  report,  we  investigate  the  effect  of  clouds 
and  precipitation  on  the  Nimbus  VI  ESMR  channel  and  the  use  of  ESMR  verti¬ 
cal  and  horizontal  components  for  the  inference  of  the  cloud  thickness  and 
the  associated  liquid  water  content  over  ocean  surfaces.  A  microwave 
radiative  transfer  program,  including  polarization,  has  been  developed  and 
comparisons  have  been  made  between  the  computed  and  observed  brightness 
temperatures  utilizing  the  available  radiosonde  and  surface  observations. 

We  then  carry  out  the  theoretical  sensitivity  analysis  and  develop  a  simple 
empirical-theoretical  equation  to  retrieve  the  cl.'ud  thickness  ar,d  liquid 
water  content. 

In  Section  3  of  this  report,  we  develop  and  derive  the  necessary 
parameterization  equations  for  infrared  and  microwave  radiative  transfer  in 
conjunction  with  the  inference  of  cloud  parameters  from  satellites.  The 
retrieval  programs  have  been  developed  for  the  determination  of  the  high 
cloud  top  height,  high  cloud  thickness,  low  cloud  top  height,  and  cloud 
liquid  water  content  successively  utilizing  selected  HIRS  and 
SCAMS  channels.  Hypothetical  error  analyses  have  been  performed,  and  it 
is  found  that  the  retrieval  technique  appears  to  be  theoretically  rigorous 
and  practically  feasible.  We  then  apply  the  cloud  sounding  technique  to 
the  Nimbus  VI  HIRS  and  SCAMS  data  for  a  number  of  carefully  selected  cases 
associated  with  summertime  convective  cloud  systems  and  wintertime  large 
scale  synoptic  cyclones  described  in  a  previous  report  by  Astling  and  Liou 
(1979)  and  a  master's  thesis  (Hall,  1979).  We  show  that  the  satellite 
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retrieved  cloud  parameters  from  a  combination  of  selected  infrared  and 
microwave  channels  appear  to  be  in  qualitative  agreement  with  the  avail¬ 
able  synoptic,  radiosonde,  surface,  and  radar  observations.  Further 
applications  of  the  cloud  sounding  technique  developed  in  this  study  for 
the  mapping  of  three-dimensional  cloud  parameters  over  the  globe  seem 
feasible. 


SECTION  2 


INVESTIGATION  OF  CLOUD  AND  PRECIPITATION  EFFECTS 
ON  THE  NIMBUS  6  ESMR  CHANNEL 

2.1  Characteristics  of  the_ESMR  Channel 

The  Nimbus  6  Electronically  Scanning  Microwave  Radiometer  (ESMR) 
receives  the  thermal  microwave  radiation  upwelling  from  the  earth’s 
surface  and  atmosphere  in  a  250  MHz  band  centered  at  37  GHz  (0.81  cm). 

The  antenna  beam  is  scanned  electrically  through  71  discrete  positions 
on  each  sweep.  Each  sweep  takes  5.3  seconds.  The  nominal  resolution 
is  20  km  crosstrack  and  45  km  downtrack.  Brightness  temperatures  are 
measured  at  each  scan  position.  The  data  are  intended  for  the  mapping 
of  the  cloud  liquid  water  content,  the  distribution  and  variation  of 
sea  ice  cover  and  snow  cover,  and  the  characteristics  of  land  surfaces. 

The  Nimbus  6  ESMR  measures  both  horizontal  and  vertical 
polarization  components  by  using  two  separate  radiometer  channels.  The 
antenna  beam  of  the  ESMR  scans  ahead  of  the  spacecraft  along  a  conical 
surface  with  a  constant  angle  of  45"  with  respect  to  the  antenna  axis 
(Fig.  2.1).  The  beam  scans  in  azimuth  '35°  about  the  forward  direction 
in  71  steps.  In  azimuth,  the  angular  resolution  varies  from  0.95° 
when  the  beam  position  is  straight  ahead  to  1.17'’  at  the  35°  scan 
extremes.  In  elevation  it  varies  From  1"  at  the  straight  forward  posi¬ 
tion  to  0.84’  at  the  extremes  of  the  scan.  The  data  are  calibrated 
using  warm  (instrument  ambient)  and  cold  (cosmic  background)  inputs  to  the 
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radiometer.  The  weight,  inn  function  for  the  37  GHz  channel  with  a  surface 
emissivity  of  0.9  is  shown  in  fig.  2.2.  At  this  frequency,  atmospheric 
absorption  is  mainly  due  to  water  vapor. 


2.2  Microwave  Radiative  Transfer  I ncl udi ng_  Pol ari zati_on 

2.2.  1  Basic  equations.  Consider  a  plane-parallel  cloud  layer  consisting  of 
absorbinn  gases  and  particulates  in  local  thermodynamic  equilibrium  and  assune 
that  thermal  radiation  from  the  earth  atmosphere  is  independent  of  the 
azimuthal  angle  We  consider,  separately,  the  radiance  components 
I  and  1  in  the  directions  parallel  and  perpendicular,  respectively, 
to  the  meridian  plane  containing  the  incident  and  emergent  beams.  The 
equations  of  radiative  transfer  governing  the  radiances  I„  and  I  can 
be  written  in  the  form 


d  I .  ( 

d,'  * 


[d  (  +  J  (,,1  )]  -  (WV)B,  [T(0], 


.(.  ,.;u)  -  'V  [J,  ( 


,r 


Jr,r( 1 ;")] 


(1-v  )B  [T(-  )], 

V 

(2.1) 


where  is  the  cosine  of  the  emergent  angle  o  which  is  with  respect  to 
the  zenith,  1  is  the  optical  depth  for  cloud  particles  and  gases  within 
the  cloud,  T  is  the  cloud  temperature,  and  the  single  scattering  albedo 
i  s 


s,- 


o  *  - 

+  K  *  k  "i  ■ 


a  ,v 


where  :  and  denote  the  volume  scattering  and  absorption  cross 

s ,  a ,  a  p 

sections,  respectively,  for  ..loud  particles  at  frequency  v,  k  the 
gaseous  absorption  coefficient,  ,  the  gaseous  density,  B  denotes  the 
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Planck  function,  and  J  ,  (1,11)  is  the  scattering  source  term  which 

P  5  P 

transforms  radiation  polarized  in  the  direction  p'  to  the  direction  p. 
Chandrasekhar  (1946)  showed  that  if  the  phase  matrix  was  determined  by 
Rayleigh  scattering,  then  the  source  terms  are 


|  I„  (v.,-i;ij,)[2(1-u*2)  +  p2(3p'2  -  2)]dp '  , 
-1 


Jr,'.'  >  ~  8  ^  Ir(o,i  ;u‘ )d|i ' 

Jv  >r('M  ;p)  =  ^  |  (v,t  ;p  ’  )ii'  2dp' 


(2.2) 


Jr,r("’i;i,) 


^  f  Ir(v»r;p')dp' 


Inserting  Eg.  (2.2)  into  Eq.  (2.1),  we  find 


dl„ (v,i ;p) 


3d>  1 


‘  =  Iv  (v, .  ;ij)  -  -8--  f/i  Ij,  (v,x;p’  )[2(1-p'2)  +  p2(3p'2 


o  1 

+  i'  ;  I r ( v , t ; i . '  )dp'l  -  (l-mv)Bv[T(T)]  , 


dl  (v,i;p)  3d)  1  o  i 

. —  -j- -  =  lr(v,t;p)  -  -g—  [  [  ;p’  )p‘  dp'  +  f  ^  lr(v,x;p' 


-  0-mv)Bv[T(7)]  • 


(2.3) 


The  solutions  of  Eq.  (2.3)  and  the  numerical  techniques  for  obtaining 
the  spatial  distribution  of  radiant  intensity  as  a  function  of  optical 
depth  by  using  a  discrete-ordinate  method  will  be  comprehensively 


2)]du ' 


)du'] 
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described  in  this  section.  For  simplicity,  we  let  I?(v,t;ii.)  "  I„  .  and 
Ir(v,r;p^)  =  1^.  Upon  replacing  the  integrals  by  the  appropriate  Gauss's 
quadrature  in  the  n^  order  approximation,  we  obtain  the  following  systems 
of  4n  linear  equations  in  the  forms 


dl 


s  orr  =  'n  -  t  {2]  +  “in  aj(3ur2)I!j 


3  « 


*  y  a,  I  01-(l-w  )B  [T(i)]  , 

j  J 
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dVi 


3(o 


i  di  Vi  “  '8  V  aj  Vj  +  4  aj"jV,j) 


(l-o»v)Bv[T(.)] 


(2.4) 

where  the  index  i  and  j  are  from  -n  to  n(n^O),  a.  represent  the  quadra- 

J 

ture  weights,  and  p..  denotes  the  emergent  angle  at  the  discrete  stream 
i.  The  complete  solutions  for  the  above  differential  equations  assuming 
an  isothermal  cloud  layer  may  be  obtained  by  seeking  a  general  solution 
for  the  homogeneous  part  plus  adding  a  particular  solution. 

2.2.2  General  solution.  The  homogeneous  part  of  Eq.  (2.4)  may  be 
written  as 


dh± 

11  i  di 


3m 


V,  -  -r  >Vj  * 


di  .  3m  9 

”i  -df  =  Vi  -  V  Vj  +  M  V.j) 


(2-5) 


In  order  to  obtain  general  sol utions  for  Eq.  (2.4),  we  let 
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(2.6) 
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Ixi  =  9i  e'kr,  lri  *  h.  e*kl  .  i  =  *1,  '2,  •••,  »n. 


where  g.'s,  h.'s  and  k  are  certain  constants.  Substituting  the 


expressions  for  Ij .  and  I  ^  into  Eq.  (2.5),  we  obtain 


3m  9  9 

(1+ijik)9i  =  r  f2^j(1'|Jj)9j +  ijitlaj(3,ij-2)9j  +  y.aj  hji}  ’ 


(Ul,,k)hi  =  r  4  ^ajhj>- 


(2.7) 


Eq.  (2.7)  implies  that  g.  and  h-  may  be  expressed  in  the  forms 


apt  +  P 


*1  ’  1‘TV;i  •  hi  *  ifsi  ■  1  '  O.  i2.  •••■  in. 


(2.8) 


where  <x,  [5,  and  i  are  certain  constants,  independent  of  i.  Inserting 


Eq.  (2.8)  into  Eq.  (2.7)  we  have 
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if  we  define  the  formula 
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then  Eq.  (2.9)  becomes 


up-  +  6  =  |  f2[u(D2-D4)  +  h{D0-D2)]  +  ,i^[u(3D4-2D2)  +  B(3D2-2Dq) 

+  yD0]  )  , 

y  =  |  (aD4  +  BD2  +  yDq)  .  (2.11) 

Since  Eq.  (2.10)  is  valid  for  all  i,  we  must  have 

|u  =  a(304-2D2)  +  B(3D2-2Dq)  +  /Dq  ,  (2.12) 


=  a(D2-D4)  +  P(DQ-D2)  ,  (2.13) 

|y  =  (xD4  +  BD2  +  y00  .  (2.14) 


Eqs.  (2.12)  -  (2.14)  represent  a  system  of  homogeneous  linear  equations 
for  a,  3,  and  /.  In  order  to  have  non-trivial  solutions  the  determinant 
of  this  system  must  vanish.  Hence 
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can  transform  the  determinant  on  the  left-hand  side  successively  into 


Hence 
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(2.16) 


(VD2-3>2  -  2<W!)(D0 -2)  =  0  • 


(2.17) 


Eq.  (2.17)  is  of  order  4n  in  k  and  admits,  in  general,  4n  distinct  non¬ 
vanishing  roots  which  must  occur  in  pairs  and 


+'t 


-  1,  2,« 


2n 


(2.18) 


From  the  relations  in  Eqs.  (2.12)  -  (2.14)  we  can  express  B  and  y  in 
terms  of  Dq,  D^,  0^,  and  <t  to  get 
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-<V°4> 

h  -  - u, 

WS 


Y  = 


P2~f3°4~D0D4 

(d0-d2-|)(d0-3) 


a  . 


Thus,  from  Eqs.  (2.6),  (2.8),  and  (2.19),  we  obtain 


4  -  (D  -D4)/(D0-02-|)  -k,T 

ijf.i  =  const‘ - T+JT.k -  e 


i  ■» 


Ir^  =  const. 


D2+'3'°4'D0D4 


-k  i 

a 


( ^ +"  i  ky ) ( Dq" ^2" 3 ) ( P0~ 3 ) 


i  =  4 1 ,4  2 , * •  •  ,fn,  <*=• 1 , 1 2 •  ,'2n. 


The  general  solution  can  now  be  written  as 


,  ?"  „  4  -  (d2-Q4)/(00-P;-3)  -kat 

(Li  a  1  i .  k  e 


c*  =  -2n 


i  i 


2n 

I  .  =  ,  M 

ri  L 


D2+3°4'D0D4 


-k  -i 

a 


a=-2n  4  (l+Miki)(D0-D2-3-)(D0-|) 


(2.19) 


(2.20) 


(2.21) 


where  M  (a=-2n,. . .2n)  are  4n  constants  of  integration. 

2.2.3  Particular  solution.  To  complete  the  solutions  for  the  system 
of  equations  denoted  in  Eq.  (2.4),  a  particular  solution  is  required. 

It  can  be  obtained  by  letting 
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Inserting  Eg.  (2.22)  into  Eg.  (2.4)  and  performing  integration  we  have 


o  O  O 

0  "  Cv.  -  8~  (!  c"  "  ^ii  +  2llicr)  '  0-«v)Bv[T(t)]. 


3o'  « 

0  =  C  -  vV  U  +  2C  )  -  (1-w  )B  [T(t  ) ]  . 
r  8  3'  .  r  \>'  v 


Eg.  (2.23)  can  be  rewritten  as 


*  4  vf»«f  '  4  4“vCr  *  0-%)Bv[T(O], 


f  C,  *  (l-y\,)Cr  -  (l-J.v)Bv[T(>n 


From  Eg.  (2.24)  we  note 


<  i  3,3  2 » r 

(1  ‘  4  V  4*V'i)Cc. 


(jh  p2  +  i  -  )c  =  o 

v  4  v  i  4  v  r 


It  follows  that 


C*  =  Cr 


Substitution  of  Eg.  (2.25)  into  Eg.  (2.24)  leads  to 


(1-%)C„_  =  (l-r„v)Bv[T(i)]  . 


(2.23) 


(2.24) 


(2.25) 


Thus,  the  particular  solution  is  simply 


(2.26) 


«■  *  i 


=  I 


i  =  Cv[T(t)]. 


2.2.4  Complete  solution.  Combining  the  solutions  given  oy  Eqs. 
(2.21)  and  (2.26),  the  complete  solutions  for  tq.  (2.4)  may  be  expressed 
in  the  form 


I 


£ ,  i 


2n  -k  i 

I  M  *  U.  )  e  u  +  B  [T( i  ) ] , 
=-2n  u  1  v 


Ir  i  =  f  M  4  (V])  e  k'1'  +  B  [T(  i )] ,  (2'?7) 

r>1  a=-2  r. 

where  ij<r  ,  and  ;!>t  represent  the  eigenfunctions,  k  the  eigenvalues  of 
the  differential  equations,  and  a  set  of  constants  of  proportionality 
to  be  determined  from  the  radiation  boundary  conditions  above  and  below 
the  cloud  layer. 

In  the  frequency  domain,  the  Planck  function  is  given  by 


Bv[T(t)]  =  2hv3/[c2(ehv/KT-l)]  , 


(2.28) 


where  h  is  the  Planck  constant,  K  the  Boltzmann  constant,  and  c  the 
velocity  of  light.  In  the  microwave  region,  the  Rayleigh-Jeans  approxi¬ 
mation  for  the  Planck  function  may  be  utilized  and  we  find  (hv/KT-+0) 

BV(T)  =  (2v2K/c2)T  .  (2.29) 

Eq.  (2.29)  states  that  the  Planck  radiance  is  linearly  proportional  to 
the  temperature.  Furthermore,  radiometers  which  measure  the  thermal 
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emission  are  usually  calibrated  with  sources  at  a  certain  reference 
temperature.  Thus,  we  may  define  an  equivalent  brightness  temperature 


Tg  such  that 

Iv  =  (2o2K/c2)Tb(v) 


(2.30) 


Using  Eqs.  (2.29)  and  (2.30),  Eq.  (2-3)  may  be  written  in  terms  of  the 
brightness  temperature  as  follows 


32i  1  , 

=  Tb,(v,i  ;,:)  -  -£■  TBp.(v,r  ;p  '  )  [2(  1 -u  ’  2 ) 

+  a 2 ( 3r 1 2~ 2 ) ]di . 1  +  „2  /  TBr( •..,Tip,)du,l  -  (1-wv)T(t), 


dTBr('.Nr 


+ 


TBr(v,T,p)  -  -g“[  [  TBfi(v,-' ;p')tJ  dp 

1 

TBr(v, .  ;u*  Jdj.’  ]  -  (1-;>v)T(t)  . 


(2.31) 


The  solutions  in  terms  of  the  temperature  field  can  be  derived  in  a 
manner  similar  to  the  procedures  outlined  previously  and  may  be 
written  in  the  form 


2n 

=  V 
/. 

u=-2n 


2n  -k  i 

V  M,'t  d>  (m  •  )e  1  +  T  , 
< =-2n  *  1  c 


(2.32) 
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where  we  note  that  the  constants  of  proportionality  M1  differ  from 
those  given  in  Eg.  (2.26)  and  are  to  be  evaluated  from  the  radiation 
boundary  conditions  in  terms  of  the  temperature  field.  Note  that  T  , 
the  cloud  temperature,  is  independent  of  frequency. 

2.2.5  Boundary  condition.  In  this  study,  it  is  assumed  that  the 
cloud  layer  is  isothermal  and  homogeneous  with  respect  to  the  cloud 
composition.  At  the  cloud  top,  the  downward  brightness  temperature  is 
equal  to  the  brightness  contributions  from  every  point  in  the  atmosphere 
above  the  cloud  top.  Thus,  we  write 


T  (v,z  i-^)  =  /  T ( z ) dV’v ( z , z t ; -|1  i ) ,  p=£,r,  (2.33) 

z=zt 

where  zt  is  the  height  of  the  cloud  top  and  -pi  indicates  that  radiation 
transfer  is  downward.  Here  we  change  the  i -coordinates  to  the  height 
coordinates  for  the  convenience  of  discussion.  The  transmittance  in  this 
equation  is  given  by 


Ty(zt,z-,-.ii)  =  T’v ( z , z t ; -v  j)  =  exp[---;  j  ky(z'  )e(z'  )dz'  ] 


(2.34) 


1  z, 


It  is  convenient  to  express  the  transmittance  in  reference  to  the  top 
of  the  atmosphere  in  the  form 


'VZt’z;iii}  =  7’v(zt’u>;ili)/7'v(z’“;‘Ji) 

=  ( 2 1 » M  i  ) / T’y  ( Z  ,P  .j  )  ,  Zt>Z  . 


(2.35) 


At  the  bottom  of  the  cloud,  two  conditions  are  considered.  The  first 
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one  is  for  a  nonpol arizing  surface  (such  as  land)  while  the  other  is 
for  polarizing  surface  (such  as  calm  water). 

for  a  nonpolar iz ing  surface,  the  radiation  incident  on  the 
cloud  base  consists  of  three  components;  (a)  the  surface  contribution, 
(b)  the  direct  atmospheric  contribution  from  below  the  cloud,  and  (c) 
the  reflected  atmospheric  contribution  from  below  the  cloud.  The 
lower  boundary  condition  can  then  be  expressed  as 


TBp('sZb;+,'i)  =  ,oVv(0’zb;+1,i)  +  I  T(z)d7’,)(z,zb;+pi) 


( 1  )y'.l(0,zb i+»‘i )  f  T(z)drv(0,z;-P1) ,  p=£,r 


(2.36) 


where  zfa  is  the  height  of  the  cloud  base. 

For  a  polarizing  surface,  we  have 

zb 

TBP(j’zb;+;,i)  =  '  wP^i)Ts‘"v(0’zb;+ui)  +  l  T(z)drv(z,zb;+Pi) 


vp(i‘i )3^(°-zb;+iii )  f  T(z)dy,v(0,z;-pi),  p=i.,r  .  (2.37) 


Each  term  on  the  right-hand  side  of  Eq.  (2.37)  has  the  same  meaning 
a^  that  in  Eq.  (2.36).  Note  that  the  emissivity  varies  with  the 
discrete  stream  employed  and  with  parallel  and  vertical  polarization 
components. 

2.2.6  Solution  matrix.  Upon  applying  the  boundary  conditions, 
Eg.  (2.32)  becomes 
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(2.38) 


'  M'<  v(t  ^ ^  =  ^  '  Fc  ’ 


=  V^t;^)  -  Tc 

u 


S  t,  <“i>  ‘  TGv^’’zb;i‘i'  -  Tc  • 

U 

l  "a  <‘1>  '  TBr(',’zb:"i>  -  Tc  • 


where 


ru(,'i)  =  e 


n  (g.-)  =  <fr  ( n  - )  e 

sX  I  vt  I 


-k  r* 

a 


-k  t* 

a 


(2.39) 

(2.40) 
(2-41) 


(2.42) 

(2.43) 


and  i*  is  the  total  optical  thickness  of  the  cloud. 

Eqs.  (2.38)  -  (2.41)  represent  a  system  of  4n  linear  equations 
from  which  the  M  coefficients  may  be  determined  by  using  standard 
matrix  inversion  techniques.  In  terms  of  matrix  operation  we  may 
express 

I  *  £  .  (2.44) 

The  coefficient  matrix  in  Eq.  (2.44)  is 
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M 

M 


-?n 

-?n+\ 


M 

-n 


rt  = 


-i 


(2.45) 


M 


n 


M 

M 


2n-l 

2n 


The  matrix  denoting  the  contributions  due  to  cloud  emission  and  upward 
and  downward  radiances  reaching  the  cloud  base  and  top,  respecti vely, 
may  be  wri tten  as 


TBp(v,zt;'"n)  '  Tc 

TB'.'(v,zt;'"l  ^  "  Tc 
TBr^v,zt;"vln)  '  Tc 

TBr(v’zt',",il^  "  Tc 
TBi'(',,zb;pl )  ‘  Tc 

TBv  ^ v,zb’|ln^  "  Tc 
TBr(v’Zb;h}  •  Tc 

TBr(v,zb;"n)  '  Tc 


(2.46) 


and  finally,  the  4n  by  4n  matrix 
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...  *_-,(-!■,) 

'  'Pr/ _‘‘n^ 

^ - 2n (  ’ 1 1 ^ 

•  -  •  'I'.-j  (-1^  ) 

'!’l(-"l) 

W'"!* 

<!,-2n|-11n) 

•••  *-l<-»*n> 

*l|-»‘n)  ••• 

1'2n^'|,n^ 

<*’-2n^~,'l  ^ 

't>_1(-p1) 

<I'1(-I'1)  •  •• 

r-2n(lJl> 

...  t_, (>.,) 

r'l{v'l) 

r2n<‘Jl> 

^-2n('Jn) 

...  >J,„) 

^(hn)  ... 

r-2nK> 

n-2n(wl) 

•  • •  n_1 (p1  ) 

n^(||.|)  ••• 

r,2n(l‘l) 

n-2n^n} 

...  rLl(pn) 

r|l(pn)  ••• 

f|2n^lln^ 

(2.47) 


Having  the  above  information,  we  may  now  carry  out  the  computations 
for  the  evaluation  of  NTs  simultaneously  by  means  of  a  matrix  inver¬ 
sion  method.  Consequently,  values  of  NTs  may  be  inserted  into  Eq. 
(2.32)  to  obtain  the  brightness  temperature  distribution. 

The  microwave  transfer  program  for  cloud  layers  will  give 
the  upward  brightness  temperature  at  the  cloud  top,  i.e.,  TD(v,z.;p  ). 

D  Li 

Thus,  the  upward  brightness  temperature  at  the  satellite  point  in 
overcast  cloudy  conditions  may  then  be  written  as 


Vv’°°;pi^  =  TB^v,ztu'i^7’v^zt’">’lii )  +  /  T(z)d7’v(z,"o;ii  ■ )  ,  (2-48) 

zt 

where  the  first  term  in  the  right-hand  side  of  Eq.  (2.48)  represents 
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the  contribution  from  the  cloud  top  brightness  which  is  being 
attenuated  to  the  top  of  the  atmosphere,  and  the  remaining  term  is 
the  atmospheric  contribution  above  the  cloud  top.  For  applications 
to  the  37  GHz  ESMR  channel,  computations  are  carried  out  for  the 
emergent  angle  of  50°. 

2.2.7  Problem  of  surface  em i ss  i v  i  ty .  A  special  problem  area  in 
the  use  of  microwave  for  atmospheric  sounding  from  a  satellite  plat¬ 
form  is  the  surface  emissivity  in  the  microwave  spectrum,  emissivity 
values  of  the  earth's  surface  vary  over  a  considerable  range,  from 
about  0.4  to  1.  Emissivity  of  the  water  surface  depends  upon  such 
variables  as  salinity,  temperature,  and  frequency.  The  brightness 
temperature  of  a  calm  sea  depends  on  the  salinity  and  temperature 
of  the  water,  the  frequency,  polarization,  and  incidence  angle. 

Over  land,  on  the  other  hand,  the  emissivity  depends  on  the  moisture 
content  of  the  soil.  Wetting  of  a  soil  surface  results  in  a  rapid 
decrease  in  the  emissivity  value.  The  surface  emissivity  appearing 
in  the  first  term  on  the  righthand  side  of  Eq.  (2.36)  has  a  signifi¬ 
cant  effect  on  the  brightness  temperature  value. 

By  assuming  an  electromagnetical ly  smooth  surface,  the 
emissivity  can  be  expressed  by 

=  1  -  r2  ,  (2.49) 

where  the  Fresnel  reflection  coefficient  for  the  air-water  interface 
is  given  by 
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(2.50) 


2  2 

(cos  ' i  -  p)  +  q 
( cos','  i  +  p)2  +  q2 

2  2  2  2 
[ (i»r  -  m^cosu.  -  p]  +  (2inrmicosOi  -  q) 

[(m2  -  m2)cosf>i  +  p]2  +  (2inrmicosOi  +  q)2 

The  subscripts  v  and  h  refer  to  the  vertically  and  horizontally 
polarized  components,  respectively.  The  angle  of  incidence  is  denoted 
as  ■!.  and  the  variables  p  and  q  are  given  by 

p2  =  fm2  -  m2  -  sin20.  +  [(m2  -  m2  -  sin2u.)2  +  4m2m2]1/2}  , 

q2  =  ~  {-(m2  -  m2  -  sin2o.)  +  [(m2  -  m2  -  sin20.)2  +  4m2m2]1/2}  , 

(2-51) 

where  mr  and  are  real  and  imaginary  parts  of  the  relative  refractive 
index  of  water  which  is  directly  related  to  the  dielectric  properties 
of  the  surface. 

2.2.8  Problems  associated  with  the  eigenvalue.  Fig.  2.3  illustrates 
a  typical  distribution  of  eigenvalues  for  Rayleigh  scattering  with  a 
single  scattering  albedo  of  0.4128.  On  the  basis  of  the  characteristic 
equation  denoted  in  Eq.  (2.17)  and  the  definition  of  D's  in  Eq.  (2.10), 
it  can  be  easily  found  that  as  k^p”\  f ( k )  ►-«-.  In  this  figure,  the 
same  intervals  between  each  p^  are  plotted  so  that  lines  across  the 
zeros  could  be  clearly  presented.  It  is  evident  that  each  interval 
contains  two  eigenvalues.  Although  the  numerical  scheme  for  the  solution 


,h‘ 


v ,  v 1 
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distribution  of  eigenv. 


of  Eq.  (2.17)  is  straightforward,  we  find  in  the  numerical  computations 
that  unstable  results  for  the  intensity  distribution  are  frequently 
generated  due  to  small  eigenvalues.  Thus,  extensive  experiments 
utilizing  Eq.  (2.17)  for  seeking  the  eigenvalue-,  corresponding  to  dif¬ 
ferent  single  scattering  albedos  have  been  carried  out  to  obtain  stable 
computations  under  various  atmospheric  conditions.  It  is  found  that 
with  the  accuracy  of  the  eigenvalues  k.  up  to  eight  decimal  points,  the 
intensities  obtained  from  the  theoretical  radiative  transfer  program  for 
the  cases  with  a  rainfall  rate  less  than  10  mm/hr  are  all  stable. 

Here  the  meaning  of  "stable"  is  used  such  that  an  increase  in  the 
accuracy  of  k^  to  about  four  decimal  points  does  not  affect  the  in¬ 
tensity  results. 

2 . 3  Comparisons  Between  Calculated  and  Observed  ESMR  Brightness 

Temperatures  for  Selected  Cases 

In  this  section  we  report  the  results  of  comparisons  between 
calculated  and  observed  brightness  temperatures  for  the  ESMR  channel. 

The  basic  criterion  for  the  selection  of  the  ESMR  37  GHz  data  for  com¬ 
parison  with  the  computed  values  is  that  the  surface  report  and  radio¬ 
sonde  data  must  be  available  under  the  satellite  pass.  Additionally,  a 
variety  of  atmospheric  conditions  and  surface  characteri sties  are  chosen 
so  that  comparisons  are  not  limited  to  specific  cases.  In  the  theoreti¬ 
cal  calculations,  a  50"  zenith  angle  was  used  in  order  to  match  the  ob¬ 
served  data  from  the  ESMR  channels.  19  cases  were  selected  at  or  near 
the  Nimbus  6  satellite  subtrack  on  August  25,  1975  and  February  20  and 
21,  1976.  The  points  were  selected  when  the  radiosonde  data  and  the 
surface  weather  report  were  available  within  30  km,  and  the  time  of  the 
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satellite  pass  was  within  six  hours  before  or  after  the  time  of  local 
sounding. 

The  Nimbus  6  has  an  equator  crossing  time  of  approximately 
1800  Z  for  the  United  States  passes.  Thus,  the  surface  report  and 
radiosonde  data  for  either  1200  Z  or  0000  Z  were  used.  We  have  collected 
11  clear  cases,  four  nonprecipitating  cloudy  cases,  and  four  precipita¬ 
ting  cloudy  cases.  Among  the  11  clear  cases,  five  were  considered  to  be 
over  polarized  surfaces  since  these  stations  were  quite  close  to  the 
ocean  surface  over  which  the  observed  ESMR  data  showed  strong  polariza¬ 
tion. 

In  order  to  have  a  reasonable  comparison  between  the  computed 
and  observed  brightness  temperatures,  the  surface  emissivity  must  be  ob¬ 
tained  objectively.  In  this  study,  we  utilize  the  scheme  proposed  by 
Waters  et  al .  (1975)  described  in  Section  3.4  to  determine  the  surface 
emissivity  which  involves  the  use  of  SCAMS  channels  1  (22.235  GHz)  and 
2  (31-4  GHz).  Since  the  water  vapor  transmittance  and  the  refractive 
index  of  water  at  37  GHz  are  very  close  to  those  at  31.4  GHz, 
we  have  assumed  that  the  surface  reflectivity  derived  for  31.4  GHz 
also  is  valid  for  37  GHz.  Consequently,  using  the  observed  ESMR 

data,  surface  emissivities  for  both  polarization  components  may  be 
estimated. 

In  Figs.  2.4  and  2.5,  the  computed  brightness  temperatures  for 
both  horizontal  and  vertical  polarization  components  are  compared  with 
the  Nimbus  6  ESMR  data.  Generally,  fairly  good  agreement  can  be  seen 
between  the  computed  results  and  observed  data.  For  those  over  polarized 
surfaces.the  computed  results  obtained  from  the  radiative  transfer 
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rig.  2.4.  Comparisons  between  calculated  and  observed  briyhtne 
temperatures  for  selected  clear  cases. 
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.5.  Same  as  Fig.  2.4,  except  Tor  cloudy  and  precipitation  cases 
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program  correctly  show  the  characteristics  of  polarization  for  the  cases 
employed.  Among  these  cases,  the  computed  horizontal  and  vortical 
brightness  temperatures  for  Miami,  and  Key  West,  Florida,  respectively, 
differ  significantly  (on  the  order  of  101'  K)  from  the  observed  value. 

The  reason  for  the  discrepancy  is  probably  due  to  the  sensitivity  of  the 
ESMR  data  on  the  surface  characteristics.  The  footprints  corresponding 
to  the  observed  data  chosen  cover  both  land  and  water  surfaces.  Hence, 
the  emissivities  derived  from  Fresnel's  reflection  coefficients  are  pro¬ 
bably  not  representing  the  actual  surface  characteristics  which  will  be 
greatly  altered  due  to  foam  and  roughness  along  the  coastline.  For 
those  over  unpolarized  surfaces,  comparison  between  the  computed  results 
and  observed  data  show  agreement  within  about  6°  K.  It  may  be  concluded 
that  the  radiative  transfer  program  developed  in  this  study  appears  to 
be  quite  reliable  despite  the  difficulty  involved  in  the  determination 
of  the  surface  emissivity  in  all  weather  conditions. 


2 . 4  Computational  Results 

2.4.1  Atmospheric  profiles  and  cloud  models.  Two  cl imatologica 1 
profiles  are  used  in  the  theoretical  sensitivity  study.  The  first 
profile  is  for  the  30°  N  latitude,  July  and  the  second  profile  is  for 
the  midlatitude  winter.  For  both  profiles,  temperature,  water  vapor 
concentration,  and  height  values  are  interpolated  to  40  standard  pres¬ 
sure  levels  for  computational  purposes.  The  temperature  and  humidity 
profiles  are  shown  in  Fig.  2.6. 

Based  on  the  foregoing  theoretical  analyses,  a  microwave 
radiative  transfer  computer  program  for  nonprecipitating  and  precipi¬ 
tating  cloudy  atmospheres  was  developed  for  both  horizontal  and  vertical 

polarization  components.  The  microwave  radiative  transfer  program  takes 
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into  consideration  effects  of  the  m  <it  ter  inn  and  absorpt  ion  ptopoit.it"-: 
of  cloud  particles  and  hydrometeors  in  a  nonhomo.jeneous  absorbing 
gaseous  atmosphere.  The  single  scatter i.ig  parameters  for  the  37  GHz 
channel  required  for  microwave  radiative  transfer  calculations  are  ob¬ 
tained  through  a  Mie  scattering  program  for  polydisperse  spheres.  The 
real  and  imaginary  indices  of  refraction  used  are  based  on  values  given 

by  Hollinger  (1973)  and  were  interpolated  to  the  37  GHz  frequency 

A  theoretical  drop  size  distribution  based  on  the  rainfall  rate 

given  by  Marshall -Palmer  (1948)  is  used  in  the  calculation.  Rainfall  rates 
of  1 ,  3,  5,  10,  15,  20,  25,  and  30  mm/hr  and  cloud  thicknesses  of  1 ,  2,  3, 

4,  and  5  km  are  used  in  the  sensitivity  analysis.  The  cloud  temperature  is 
assuned  to  be  the  same  as  that  of  the  surrounding  air  and  a  mean  temperature 
value  is  taken  in  conjunction  with  the  radiative  transfer  calculation.  To 
carry  out  theoretical  calculations  over  the  land  surface,  ari  emissivity  value 
of  0.96  is  used  for  each  discrete  stream,  lor  the  water  suit.-Ke.  the 
refractive  index  of  the  air-water  interface  depends  on  the  temperature 
and  salinity.  Shown  in  Table  2.1  are  emissivity  values  for  a  number  of 
surface  temperatures  and  salinities.  Since  variations  in  the  emissivity 
value  are  quite  small,  we  have  assumed  a  temperature  of  273u  K  and 
35  PPT  (percentage  per  thousand)  salinity  in  the  theoretical  computations. 
Fig.  2.7  shows  the  degree  of  polarization  over  the  water  surface  as  a 
function  of  the  incident  (or  emergent)  angle.  It  is  seen  that  there  are 
differences  in  values  between  these  cases  but  the  polarization  patterns 
are  all  quite  similar.  Over  the  water  surface,  the  emissivity  value  for 
each  discrete  stream  has  been  calculated  by  using  Eqs.  (2.49)  -  (2.51). 

2.4.2  Unpolarized  surface  (over  land).  Brightness  temperatures  for 
both  horizontal  and  vertical  polarization  components  over  the  land  sur¬ 
face  as  functions  of  the  discrete  zenith  angle,  and  cloud  thickness  for 


1 

( 
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Salinity  35  PPT 


INCIDENT  ANGLE  (deg) 


Defjffv  oi  I  inpdr  |>o I nr  i/at.  ion  over  the  water  surface  as  a 
function  of  the  incident  (emergent)  angle. 


Table  2.1.  Emissivi ties  over  water  surface  at  50"  incidence 
angle. 


Temperature 

(°C) 

Salinity 

PPT 

Refractive  Index 

"V  mi 

Emiss 

'  h 

ivity 

'  V 

Degree  of 
Polarization 
(%) 

0 

35 

3.8862 

2.3838 

0.3796 

0.6811 

28.90 

10 

35 

4.3893 

2.6158 

0.3495 

0.6502 

30.08 

20 

35 

4.9625 

2.7670 

0.3256 

0.6181 

30.99 

30 

35 

5.4878 

2.7849 

0.3105 

0.5972 

31.58 

0 

0 

3.8572 

2.3189 

0.3853 

0.6951 

28.67 

0 

10 

3.8748 

2.3810 

0.3802 

0.6889 

28.87 

0 

20 

3.8848 

2.3844 

0.3796 

0.6882 

28.90 

0 

30 

3.8875 

2.3848 

0.3795 

0.6398 

25.54 

rainfall  rates  are  of  1  and  10  nm/hr  displayed  in  Fig.  2.8  and  2.9. 

The  temperature  and  humidity  profiles  for  30  N  latitude,  July  and  an 
emissivity  value  of  0.96  for  each  discrete-stream  were  used  in  the 
theoretical  calculations.  Brightness  temperatures  emerging  from  the 
rain  clouds  over  land  surface  show  weak  polarization  as  evident  in  these 
figures.  Moreover,  it  is  seen  that  the  brightness  temperature  decreases 
as  the  cloud  thickness  and  rainfall  rate  increase.  Clearly,  over  land 
surfaces,  the  total  liquid  water  content  in  the  atmosphere  reduces  the 
upwelling  brightness  temperature. 

Figure  2.10  depicts  variations  of  the  brightness  temperature  at 
50°  incidence  angle  emerging  from  the  top  of  the  atmosphere  as  functions 
of  the  cloud  thickness  and  rainfall  rate.  Brightness  temperatures  de¬ 
crease  when  the  cloud  thickness  increases  as  pointed  out  previously. 

When  the  rainfall  rate  is  less  than  about  5  mm/hr,  brightness  tempera¬ 
tures  decrease  drastically  as  rainfall  rates  increase  for  all  five  cloud 


thicknesses.  The  decrease  is  much  more  pronounced  for  thicker  clouds. 
We  see  that  brightness  temperatures  do  not  change  significantly  when 
the  rainfall  rate  is  greater  than  5  mm/hr.  The  reason  that  the  bright¬ 
ness  temperature  decreases  when  the  rainfall  rate  increases  is  due  to 
the  progressively  stronger  backscattering  by  large  raindrops.  Note 
that  in  reality  thin  clouds  normally  do  not  produce  high  rainfall  rates 
2.4.3  Polarized  surface  (over  calm  water).  In  this  subsection,  we 
report  the  resulting  calculations  which  were  carried  out  for  polarized 

surfaces.  Let  the  horizontal  and  vertical  brightness  temperatures  be 
H  V 

Tg  and  Tg,  respectively  and  define  the  degree  of  linear  polarization 
as  follows: 


P(%) 


T^  - 

'b  'b 

TV  +  TH 

'b  'b 


x  100% 


(2.53) 


The  polarized  emissivities  used  in  the  transfer  calculation 
were  derived  from  Eqs.  (2.49)  -  (2.51)  using  the  index  of  refraction 
corresponding  to  a  temperature  of  273°  K  and  a  salinity  of  35  PPT. 
Emissivity  values  as  functions  of  the  incident  (or  emergent)  angle  for 
both  horizontal  and  vertical  components  are  shown  in  Fig.  2.11.  The  sur¬ 
face  polarization  varies  from  0%  at  the  0°  incident  angle  to  90%  at  the 
1  imb. 

In  Figs.  2.12  and  2.13  are  shown  both  vertical  and  horizontal 
polarization  components  over  the  water  surface  as  functions  of  the  in¬ 
cidence  angle  using  a  30°  N,  July,  profile  for  1  and  10  mm/hr  rainfall 
rate.  The  brightness  temperature  patterns  over  polarized  surfaces  show 
significant  differences  from  those  over  unpolarized  surfaces. 
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SURFACE  EMISSIVITY 


July  /  mm/hr 


Brightness  temperatures  for  both  horizontal  and  vertical 
components  over  ocean  surface  as  t'uer  ^ior.c  of  +he  °mercen 


Significant  differences  between  horizontal  and  vertical  components  at 
the  top  of  the  atmosphere  can  bo  seen  for  low  rainfall  rates  and  thin 
clouds  over  the  water  surface.  Moreover,  increasing  the  rainfall  rate 
and  cloud  thickness,  the  brightness  temperature  increases  for  emergent 
angles  less  than  about  60°  owing  to  the  larger  cloud  emissivity  as  com¬ 
pared  with  the  surface  emissivity.  The  polarization  characteristics  of 

the  underlying  surface  are  found  to  have  little  effect  on  the  emergent 

brightness  temperature  for  high  rainfall  rates.  For  a  rainfall  rate  of 

1  mm/hr,  each  cloud  model  produces  quite  different  brightness  tempera¬ 
ture  distribution  between  two  polarization  components,  especially  for 
the  1  km  cloud  case.  High  degree  of  polarization  is  seen  to  be  located 
between  40°  -  70°  incidence  angles  (Fig.  2.14).  As  the  cloud  becomes 
thicker,  polarization  is  found  to  be  greatly  reduced.  In  Figs.  2.15 
and  2.16,  brightness  temperatures  emergent  from  a  midlatitude  winter 
atmosphere  are  illustrated.  In  general,  the  brightness  temperature 
patterns  for  vertical  and  horizontal  components  are  found  to  resemble 
those  from  a  30°  N,  July,  profile.  As  can  be  seen  in  Fig.  2.14,  the 
emergent  polarization  values  for  two  different  profiles  are  rather  close. 
Apparently,  effect  of  the  atmospheric  temperature  profile  on  the  emer¬ 
gent  polarization  pattern  appears  to  be  minimum. 

Figures  2.17  and  2.18  show  the  variation  of  the  brightness 
temperature  at  the  top  of  the  atmosphere  as  a  function  of  the  rainfall 
rate  for  the  30°  N,  July,  and  midlatitude  winter  profiles,  respectively. 
Both  figures  reveal  quite  different  results  from  those  over  the  un¬ 
polarized  surface.  For  a  given  cloud  thickness,  the  brightness  tempera¬ 
tures  for  both  polarization  components  increase  with  increasing  rainfall 


41 


Degree  of  polarization  as  a  function  cf  emergent  angle  for  five  cloud 
thicknesses  having  a  rainfall  '"ate  of  1  mm/hr. 
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Same  ar  Fig.  2.12,  except  using  a  Mid-latitude  winter  profile. 


BRIGHTNESS  TEMPERATURE  (•  R ) 


rate.  The  horizontal  component  increases  more  rapidly  as  the  rainfall 
rate  increases.  For  a  1  km  thick  cloud,  the  difference  between  rainfall 
rates  of  1  inn/hr  and  3  mm/hr  is  30"  K,  but  it  is  only  3‘  K  between  25 
min/hr  and  30  mm/hr.  For  the  vertical  component  we  see  a  5°  K  difference 
between  1  mm/hr  and  3  mm/hr  rainfall  rates.  Comparing  the  brightness 
temperatures  for  cloud  thicknesses  1  and  5  km,  it  is  seen  that  the  dif¬ 
ference  between  1  mm/hr  and  3  mm/hr  is  20°  K  for  a  1  km  thick  cloud, 
whereas  it  is  44°  K  for  a  5  km  thick  cloud  for  the  horizontal  component. 
Only  5°  K  and  9°  K  differences,  respectively,  are  observed  for  the  verti¬ 
cal  component.  Fig.  2.19  illustrates  the  degree  of  linear  polarization 
as  a  function  of  the  rainfall  rate  for  the  two  atmospheric  profiles  at 
an  incident  angle  of  50°.  The  polarization  patterns  are  fairly  similar 
for  the  two  different  profiles.  Degree  of  polarization  at  this  angle  is 
seen  to  be  reduced  as  the  rainfall  rate  and  cloud  thickness  increase. 

That  is  that  the  total  liguid  water  content  within  the  atmosphere  has  a 
profound  effect  on  the  emergent  polarization  configuration  at  the  37  GHz 
frequency.  In  a  cloudy  or  precipitating  atmosphere,  which  is  over  a 
water  surface,  the  emergent  polarization  varies  from  0 1  to  about  20%. 
Because  of  these  variations,  it  seems  possible  to  derive  information  of 
the  cloud  thickness  and  rainfall  rate  (or  liquid  water  content)  over 
water  surfaces  from  vertical  and  horizontal  polarization  components.  The 
relationship  between  the  rainfall  rate  and  liquid  water  content  is 
illustrated  in  Table  3.1. 

2 . 5  Liquid  Mater  Content  and  Cloud  Thickness  Determination  from 

ESMR  37  GHz  Channel 

On  the  basis  of  the  theoretical  analysis,  it  appears  quite 
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clear  that  the  cloud  thickness  and  rainfall  rate  (or  liquid  water  content) 
in  the  atmosphere  must  be  related  to  the  emergent  vertical  and  horizontal 
brightness  temperatures.  The  functional  relationship  between  the  bright¬ 
ness  temperatures  and  cloud  parameters  may  be  rather  involved  and  compli¬ 
cated.  However,  in  an  attempt  to  understand  how  well  the  37  GHz  polari¬ 
zation  data  could  be  used  for  the  quantitative  inference  of  cloud  and 
precipitation  properties  over  ocean  surfaces,  it  is  assumed  that 


Tg  =  aQ  +  a-|AZ  +  a2AW 
Tg  =  bQ  +  b-jAZ  +  b2AW 


(2.53) 


where  ai  and  bi  (i=0,l,2)  are  certain  relevant  constants.  It  follows 
that 


AW  ^  wQ  +  +  w2tJ 

AZ  z0  +  zlTg  +  z2Tb 

where  w.  and  z.  (i=0,l,2)  are  coefficients  to  be  determined  through 
theoretical -empirical  procedures . 

A  three-dimensional  linear  regression  program  was  utilized  to 
determine  the  coefficients  w^  and  z^  denoted  in  Eqs.  (2.53)  and  (2.54) 
over  calm  water  surfaces.  A  total  of  31  atmospheric  cases  were  used  for 
30°  N  July  and  midlatitude  winter  profiles  separately  to  compute  bright¬ 
ness  temperatures  using  the  microwave  radiative  transfer  program  described 


(2-54) 
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in  the  previous  sections.  Different  1 iquid  water  content  and  cloud 
thickness  were  inserted  into  the  atmosphere  to  obtain  their  relationship 
with  the  computed  brightness  temperatures  T"  and  TZ.  The  empirically 

D  D 

derived  coefficients  based  on  the  theoretical  data  are  listed  in  Table 

2.2. 

In  order  to  examine  the  validity  of  the  empirical  formulas, 
a  number  of  theoretically  computed  brightness  temperatures  were  used 
to  retrieve  the  liquid  water  content  and  cloud  thickness.  The  true 
and  retrieved  values  are  listed  in  Table  2.3  along  with  the  brightness 
temperatures,  for  large  liquid  water  content,  the  retrieved  results 
appear  to  be  within  the  acceptable  range.  For  small  liquid  water  con¬ 
tent,  negative  values  are  obtained  through  the  simple  empirical  equations. 
In  view  of  the  extreme  difficulties  encountered  in  the  sounding  of  cloud 
parameters  from  satellites,  the  retrieved  results  derived  from  the  theo¬ 
retical-empirical  analysis  seems  encouraging. 

Thus,  we  have  selected  a  number  of  ESMR  data  over  ocean  areas 
to  test  this  technique  for  the  cloud  parameter  inference.  The  ESMR  data 
chosen  cover  mainly  the  area  20-25°  N  and  70-80°  W  for  August  25  cases 
and  20-25°  N  and  80-85°  W  for  February  20  cases  described  in  a  report 
by  Astling  and  Liou  (1979).  The  areas  under  consideration  on  these  days 
had  significant  cloud  cover.  Results  of  the  retrieved  liquid  water  con¬ 
tent  and  cloud  thickness  are  tabulated  in  Table  2.4.  In  these  exercises, 
empirical  equations  derived  from  the  30°  N,  July  and  midlatitude  winter 
profiles  were  utilized  for  August  and  February  cases,  respectively. 

Results  are  surprisingly  reasonable,  except  for  a  thin  cloud  case 
where  a  negative  liquid  water  content  is  obtained.  Flowever,  because  of 
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Table  2.2  Empirical  coefficients  for  the  determination  of 
liquid  water  content  and  cloud  thickness  over 
water  surfaces. 


Profile 

wo 

W1 

W2 

Z0 

Z1 

Z2 

30°  N,  July 

1.0428 

6.4866 

-9.8447 

1 .3747 

1 . 6896 

-5.9171 

xlO"3 

xlO"3 

xlO2 

xlO'1 

xlO"1 

Midlatitude 

1-0850 

6.9790 

-10.7565 

1 .1770 

2.3442 

-7.4291 

winter 

xlO'3 

xlO'3 

xlO2 

xlO"1 

xlO'1 

the  lack  of  surface  weather  observations  over  ocean,  it  is  not  possible 
to  verify  the  degree  of  realiability  of  these  retrieval  analyses.  Cer¬ 
tainly,  remote  sensing  techniques  developed  for  the  cloud  inference  re¬ 
quire  proper  verification  from  simultaneous  in-situ  aircraft  cloud 
physics  measurements  and  surface  observations. 
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Table  2.3.  Hypothetical  analyses  for  the  cloud  liquid 
water  content  and  thickness  retrieval;  (a) 
30°  N,  July  profile,  and  (b)  midlatitude 
winter  profile. 


(a) 


Actual 

Brightness 

Retrieved 

AW  _ 

AZ 

Temperature 

(°K) 

AW  - 

AZ 

( gni/cni  ) 

(km) 

H 

V 

(gm/cnr) 

(km) 

0.02668 

3 

163.17 

220.61 

-0.0707 

2.10 

0.02238 

1 

169.04 

223.20 

-0.0580 

1 .55 

0.06874 

2 

210.02 

233.47 

0.1066 

2.39 

0.30766 

5 

237.37 

234.74 

0.2715 

6.26 

0.34002 

4 

241  .75 

238.74 

0.2606 

4.64 

0.33043 

3 

246.16 

242.96 

0.2477 

2.88 

0.02668 

3 

(b) 

146.01 

202.29 

-0-0719 

1.65 

0.02238 

1 

149.80 

203.67 

-0.0603 

1 .51 

0-06874 

2 

193.33 

216.54 

0.1050 

2.15 

0.30766 

5 

221.73 

219.28 

0.2737 

6.78 

0.34002 

4 

226.25 

223.40 

0.2606 

4.75 

0.33043 

3 

230.63 

227.63 

0.2461 

2.66 

Table  2.4.  Retrieved  liquid  water  content  and  cloud 
thickness  for  (a)  August  25  cases,  and  (b) 
February  20  cases 

Location  Observed  Retrieved 


Lat. 

Long. 

Brightness 

Temperature 

(°K)  AW  „ 

AZ 

(°N) 

(°W) 

H 

V 

(gm/cm; 

(km) 

22.8 

80.3 

192 

223 

0-0928 

5.54 

22.8 

83.6 

180 

221 

0.0347 

4.70 

22.8 

80.5 

189 

219 

0.1127 

7.40 

22.8 

79.8 

174 

219 

0.0154 

4.87 

21.6 

79.8 

169 

217 

0.0027 

5.21 

19.5 

79.3 

172 

213 

0.0615 

8.08 

22.2 

81 .0 

155 

Tb) 

200 

0.0154 

5.46 

22.1 

84.4 

170 

208 

0.0341 

3.03 

23.8 

81.5 

169 

206 

0.0486 

4.28 

23.8 

81.9 

166 

204 

0.0492 

5.06 

43.6 

86.5 

145 

202 

-0.0759 

1.63 

SECTION  3 


EXPLORATION  OF  CLOUD  PARAMETER  INFERENCE  FROM  A 
COMBINATION  OF  INFRARED  AND  MICROWAVE  CHANNELS 

3 . 1  Par ameteri zation  of  Infrared  Rad iative  Transfe r  for 

a  Two- Layered  Cloud  System 

In  the  infrared  15  um  CO2  region,  the  upwelling  radiance  in 
a  cloud-free  atmosphere  may  be  written  as 

a- 

I -H  =  B-(Ts)r-K0)  +  /  B-[T(z)]  K-(«  ,z)dz,  (3.1) 

where  B-  denote  the  Planck  function  at  wavenumber  ">,  T  the  surface 

v  s 

temperature,  T-  the  transmittance,  and  K-  the  weighting  function. 

We  first  consider  an  overcast  single  cloud  layer  in  the 
atmosphere.  Since  the  cloud  reflectivity  is  generally  quite  small  in 
the  infrared  region  (Liou  et  al . ,  1978;  Hunt,  1973),  the  effect  of  cloud 
reflection  may  be  ignored.  Thus,  the  basic  equation  for  the  upwelling 
radiance  may  be  written  as 

I;H  "  <Az>  +  ^¥T(zc>^(~’zc> 


+  /  B-[T(z)]K-Kz)dz  ,  (3.2) 

zc 

where  z(  is  the  cloud  top  height,  7^  and  1^  represent  the  cloud  trans¬ 
missivity  and  emissivity,  respectively,  and  I-(z  )  is  the  upwelling 
radiance  reaching  the  cloud  bottom  which  is  given  by 
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(3.3) 


I,-(zc)  ■  B-(Ts)',--(zc,0)  *  /V[T(z)]K-(zc,z)dz  . 

In  Eq.  (3.2).  the  first  term  on  the  right  hand  side  represents  the 
radiance  contribution  front  the  atmosphere  below  the  cloud,  the  second 
term  is  the  contribution  due  to  cloud  emission,  and  the  final  term  is 
the  contribution  of  the  gaseous  emission  from  the  atmosphere  above  the 
cloud.  In  Eq.  (3.3),  ''’-(z  ,0)  represents  the  clear  column  transmittance 
in  reference  to  the  cloud  top  height  z  .  However,  the  clear  column 
transmittance  are  normally  available  in  reference  to  the  top  of  the 
atmosphere;  i.e.  in  the  form  of  '/’-(»', 0).  Hence,  it  is  convenient  to 
change  the  former  transmittance  to  the  latter  form.  According  to  the 
definition  of  the  spectral  transmittance,  it  is  the  ratio  of  the  output 
radiance  to  the  input  radiance.  Let  the  input  radiance  for  the 
first  layer  be  IQ  and  the  output  radiance  be  1^  ,  which  is  also  the 
input  radiance  for  the  second  layer.  If  the  output  radiance  for  the 
second  layer  is  I^,  then  the  transmittance  of  the  combined  layer,  by 
definition,  can  be  obtained  from  the  product  of  the  transmittance  of 
the  individual  layer,  i.e.,  I  1  g  =  ( Ig/I-i  )0]/ Iq)  ‘  *3as^s 

this  physical  argument,  the  transmittance  at  a  given  height  below  the 
cloud  layer  in  reference  to  the  top  of  the  atmosphere  may  be  expressed 
by 


Thus,  the  transmittance  contained  in  the  surface  term  of  Eq.  (3.3)  may 
then  be  written  as 
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(3.4) 


T~{ Zc,0)  =  T-h,0)/T-{™,Zr)  , 

and  the  weighting  function  is  now  given  by 

w-(zc,z)  i  z)  , 

Kv(zc’z)  Dz  "  r^VzT  l)z  "  r^-VzT  M0,,z)  • 

(3.5) 

For  simplicity,  K-(<»,z)  will  be  simply  denoted  as  K-(z)  in  the  follow¬ 
ing  analysis.  Upon  substituting  Eqs.  (3.4)  and  (3.5)  into  Eg.  (3.3) 
and  Eq.  (3.2),  the  upwelling  radiance  for  an  atmosphere  containing  a 
single  cloud  layer  is  now  given  by 

z 

r  r  C 

I-H  =  !T- ( Az ) { B- ( T  ) T- («> , 0 )  +  {  B-[T(z)]K-(z)dz) 

V  V  v  S  V  v  V 

(JO 

+  e£(Az)B-[T(zc)]7’-Kzc)  +  /  B-[T(z)]K-(z)dz  .  (3.6) 

zc 

We  now  consider  an  atmosphere  containing  multiple  cloud 

layers.  Following  the  foregoing  analyses  for  a  single  cloud  layer,  we 
he  1c 

let  T-  and  t-  denote  the  transmissivity  and  emissivity  for  high  and 

V  V 

low  cloud,  respectively.  Then  the  upwelling  radiance  of  an  atmosphere 
with  a  high  cloud  layer  can  be  written  in  a  form  similar  to  Eq.  (3.2) 
as  follows: 


l£(.)  -  f^tAZh)l-(zh)+ebC  B-  CT(Z „)])!•;(<»,*„> 


+  /  B-[T(z)]K-(°°,z)dz  ,  (3.8) 

zh 
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where  zh  and  Az()  are  the  high  cloud  top  height  and  thickness,  respectively. 
Assuming  that  there  is  a  low  cloud  deck  under  the  high  cloud,  then  the 
upwelling  radiance  reaching  the  high  cloud  base  is  given  by 


!-(zh)  =  n-(z,,)'/f 


(AZ£>  +  'f(A^)B-[T(zv)]}7’-(zh,Z?) 


zh 

+  f  B-[T(z)]K-(zh,z)dz  , 
z„ 


(3-9) 


where  zv  and  A z (/  are  the  low  cloud  top  height  and  thickness,  respectively, 
and  I-(z,, ),  the  upwelling  radiance  reaching  the  low  cloud  layer,  is 


I-(z,,)  "  B-(Ts)7'-(2„,0)  +  f  B-[T(z)]K-(zjl,z)dz  . 


(3.10) 


Since  the  cloud  reflection  in  the  infrared  15  urn  CO^  band  is 
normally  very  small  (see  Liou  et  al,  1978),  we  may  write 


(3.11) 

^C(Aze)  =  1  -  7f(Azf) 

Upon  substituting  Egs.  (3.10)  and  (3.11)  into  Eq.  (3.9),  we  find 

z 

Mzh)=  fBv(T5)?v(zi’°)  +  f  ^[T(z)]K-(zv:')dz}T^c  (Az£)T-(zh  ,ze ) 

+  [l-'^C(AzJ  ]B-[ T ( z . , )  ] 7'-(z h , z „ ) 
zh 

*  f  B-[T(z)]K-(zh,z)dz  .  (3.12) 
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We  then  change  the  reference  layer  of  clear  column  transmittances  and 
weighting  functions  from  high  to  low  cloud  layer  to  the  top  of  the 
atmosphere  as  described  previously  to  obtain 


Tv{h'0)  =  r-(“,°)/[j*-(».2h)y-(zh,z(l)] 


K-(zrz)  =  K-1(«,z)/[3,;:(“»zh)r-(zh,z„)]  ,  0<z<z(l, 


’  h'‘v'  h  t 


K-(zh,z)  =  K-(»,z)/r-(“’,zh)  ,  Z„  •  Z  r  z 


(3.13) 


Substituting  Eq.  (3.13)  into  Eq.  (3.12)  and  subsequently  into  Eq.  (3.8), 
we  obtain  the  upwelling  radiance  emergent  from  a  two-layer  cloudy  atmos¬ 
phere  in  the  form 


I °-H  =  T^(Mh)7^C(M9){B-(ls)T-h,0)  +  /  B-[T(z)]K-(z)dz 


v *•  '  J  v 


-  B-[T(z!i)]r-(™,zs,)}  +  T*C{ Azh)  f  /  B-rT(z)]K-(z)dz 


+  Bv£T(z*>]^<".zi>  -  B0(zh)7,v(a5’zh)}  +  BvtT(zh)]Tv(w’zh) 


+  /  B-[T(z)]K-(z)dz  • 


In  the  absence  of  high  cloud,  7^c(Azh)  =  1,  and  zh=z? , 


(3.14) 


Eq.  (3.14)  reduces  to 


p 

I^CH  -  yf  (az,);b-(ts)7'-(-,0)  +  (  B-[T(z)]K-(z)dz} 

CO 

+  [l-7'J-;C(Az(;)]B-(zfc)T-(-'>,z?)  +  /  B-[T(z)]K-(z)dz  .  (3.15) 

z 

This  is  the  upwelling  radiance  for  the  single-layer  cloudy  atmosphere. 
Likewise,  if  the  low  cloud  is  absent,  then  Eg.  (3.14)  becomes 

zh 

,?(-)  *  ¥Ts>^’0>  +  f  B-[T(z)]K-(z)dz] 

+  [l-7^c(Azh)]B-(zh)r-(«  ,zh)  +  f  B-[T(z)]K-(z)dz  .  (3.16) 

zh 

In  no  cloud  conditions,  7 -|C  ( A  z  ^ )  =  7'-c(Azv)  =  1,  and  z^  =  z^  =  0, 

Eq.  (3.14)  reduces  to  the  expression  for  the  clear  column  radiance 
denoted  in  Eq.  (3.1).  Eg.  (3.14)  is  the  fundamental  parameterized 
transfer  equation  to  be  used  for  the  inference  of  the  cloud  thickness, 
and  high  and  low  cloud  top  heights . 

3.2  Hipjh  and  Low  Cloud  Top  and  Thickness  Determination 

from  HIRS  Channels 

3.2.1  Determination  of  cirrus  cloud  top  height.  In  the  following 
analysis,  we  replace  "  by  i  to  indicate  the  channel  number  in  the  15  pm 
CO^  band.  Next  we  let  the  difference  between  the  upwelling  radiances 
for  atmospheres  containing  no  cloud  and  a  single-layer  high  cloud  layer 
be 


A  A. 

i 


i  (”)  -  iJcH 


(3.17) 
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Substituting  Eqs.  (3.1)  and  (3.16)  into  Eq.  (3.17),  we  have 


AAi  =  [l-rJc(Azh)]lBi(Ts)ri(-..,o)  -  Bi[T(zh)]ri(-,zh) 
zh 

+  /  B.[T(z)]K.(z)dz)  (3.18) 

0 


To  a  good  approximation,  the  cloud  transmissivity  may  be  considered 

to  be  independent  of  the  channel  number  in  the  15  pm  interval  (Liou, 

et  al.,  1978),  i.e.,  T^(Az)  =  TC (Az) ,  where  i  and  j  are  two  adjacent 

channels  in  the  15  pm  CO^  band.  Define  the  ratio  of  the  radiance 

differences  for  i  and  j  channels  as  follows: 

zh 


G(zh^  AA 


AA.  9i<2h>  +  l  Bi[T(z)]Ki(z)dz 


2h 


(3.19) 


J  9j(zh)  +  j  B.[T(z)]K.(z)dz 


where 


gk(zh}  =  Bk(Ts)7kKo)  ’  Bk[T(zh)]rk(a,’zh)  ’  k=Uj  ■ 

Because  of  the  ratio,  the  G  function  is  independent  of  the  cloud 
opacity,  but  depends  only  on  the  weighting  function  of  the  channel,  the 
cloud  height,  and  the  radiance  profile  from  the  surface  to  the  cloud 
top  (McCleese  and  Wilson,  1976). 

The  G  function  has  been  evaluated  based  on  Eq.  (3.19)  utilizing 
the  channel  pair  (2,3),  (3,4),  (4,5),  and  (5,6)  of  the  Nimbus  6  HIRS 
15  pm  band  for  two  climalogical  temperature  profiles.  In  Fig.  3.1, 
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solid  .Hid  dashed  1  inns  are  for  mid-latitude  summer  and  a  tropical 
summer  profile,  respectively.  It  is  seen  that  the  6  function  is  not 
strongly  dependent  on  the  atmospheric  temperature  profile  up  to  the 
level  of  the  tropopause.  Inspection  of  Fg.  (3.19)  reveals  that  the 
G  value  depends  only  on  z^  for  a  given  radiance  profile  B-[T(z)]. 

The  G  function  normally  varies  insignificantly  when  the  difference  of 
the  radiance  profiles  is  not  very  pronounced.  The  relation  between 
the  G  value  and  z^  may  he  constructed  through  a  regression  analysis. 

T rum  the  definition  of  t ho  G  function,  the  observed  and  clear  column 
radiance  values  are  required  for  the  computation.  These  values  may 
be  evaluated  from  the  observed  temperature  profile. 

In  Fig.  3.1,  curves  A,  B,  C,  and  D  are  for  the  channel  pair 
(2,3),  (3,4),  (4,5),  and  (5,6),  respectively  for  a  mid-latitude  summer 
profile.  We  find  that  the  regions  in  which  the  cloud  height  z ^  is  a 

monotonia  function  of  the  G  value  are  quite  limited  for  curves  A  and 

P>.  We  note  that  the  spread  of  the  G  value  is  from  0.03  to  0.08, 
while  z ^  ranges  from  10  to  12.8  km  for  curve  A.  This  implies  that  for 
any  G  value  outside  the  domain  from  0.03  to  0.08,  there  will  be  no 
solution  for  z,  ■  Moreover,  unless  z.  is  between  10  and  12.8  km,  there 
also  will  be  no  solution  for  it.  The  region  in  which  the  G  value  in¬ 
creases  with  height  is  much  wider  for  curves  C  and  D.  The  G  value 

ranges  from  0.02  to  0.32  when  z h  is  between  2.4  and  12  km  for  curve  C. 
The  slope  of  the  curve  C  is  less  charp  as  compared  with  curves  A  and  8.  For 
this  curve  z^  is  less  affected  by  the  change  of  the  G  value.  This 
implies  that  the  estimated  z^  value  is  more  tolerable  to  the  observed 
error-  in  the  G  value. 
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3.1.  The  relationship  between  the  b  function  (.ee  text  for  the 
definition)  and  the  hiqh  cloud  top  heiqht  for  four  channel 
pairs  us  ini)  two  atmospheric  profiles. 
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In  view  of  the  proceeding  consideration,  we  find  that  curves 
C  and  D  are  better  options  for  the  determination  of  the  high  cloud  top 
height.  Since  channel  6  senses  that  the  lower  part  of  the  atmosphere, 
the  low  cloud  has  a  significant  effect  on  the  determination  of  the  G 
value  for  the  channel  pair  (5,6).  In  addition,  the  region  for  the 
curve  D  extends  to  only  about  11  km.  It  follows  that  this  channel  pair 
is  not  ideal  for  estimating  the  high  cloud  top  height. 

for  curve  C  which  is  constructed  from  channels  4  and  5,  the 
cloud  height  z^  is  approximately  a  linear  function  of  G  between  6  and 
12  km.  The  height  range  is  approximately  the  location  of  the  cirrus 
cloud  top.  Thus,  z^  may  be  determined  by  the  following  linear  equa¬ 
tion  : 


where  the  coefficients  d^  and  d^  are  determined  through  a  regression 
analysis. 

3.2.2  Determination  of  cirrus  cloud  thickness.  The  HIRS  water 
vapor  channel  in  the  6.3  cm  H^O  band  is  intended  for  the  estimate  of 
the  water  vapor  concentration  in  the  upper  troposphere  and  the  detec¬ 
tion  of  thin  cirrus.  The  weighting  functions  for  the  HIRS  temperature 
and  water  vapor  channels  are  shown  in  Tig.  3.2  (Smith  et  al.,  1975). 
Examination  of  the  weighting  function  of  the  6.7  urn  water  vapor  channel 
(channel  10)  reveals  that  this  channel  peaks  at  about  400  mb  and  has 
very  insignificant  interference  from  the  atmosphere  below  700  mb. 

After  the  cloud  top  height  is  estimated  from  the  previous  analyses  the 
cloud  thickness  a Kc;  may  be  evaluated  from  the  scheme  described  below. 
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Tor  a  tliin  cirrus  cloud,  the  cloud  transmissivity  may  be 
approximately  expressed  as  an  exponential  function  of  the  cloud  optical 
thickness.  Since  i/  (km  is  the  total  extinction  of  radiation  due  to 

absorption  and  scattering  by  the  cloud  particles,  the  geometric 
thickness  Az(km)  is  related  to  the  optical  depth  by  the  equation 


Az  =  i /b  -  (kill) 


(3.20) 


For  the  computation  of  [•>  ,  the  cirrus  cloud  was  assumed  to 
be  composed  entirely  of  ice.  The  single  scattering  computation  was 
made  for  the  central  wavenumber  1508.29  cm  ^  in  the  water  vapor  band 
utilizing  a  theoretical  ice  cloud  model  developed  by  Liou  (1972).  The 
major  and  minor  axes  of  the  cylinder  used  were  200  and  60  pm,  respec- 

_3 

tively,  with  a  concentration  of  0.05  cm  .  The  real  and  imaginary 
parts  of  the  index  of  refraction  were  taken  from  the  measurement  by 
Schaaf  and  Williams  (1973).  The  value  for  ft  is  1.326  km  ^ . 

From  Eg.  (3.20),  the  transmissivity  of  a  thin  cirrus  cloud 
may  be  expressed  as 


"i  =  ex[,(-  l^e  Azh}  ' 


(3.21) 


Upon  substituting  Eq.  (3.21)  in  Eq.  (3.16),  the  cirrus  cloud  thickness 
may  be  derived  to  give 


-  f f’i (zh)'/’i (--,zh)  +  I  bi[T(z)]Ki(z)dz) 


/  Bi[T(z)]Ki(z)dz  -  Bi(zh)7’i(c"’Zh) 


(3.22) 


Because  of  the  strong  water  vapor  absorption  in  this  channel,  low 

clouds  below  700  mb  have  little  effect  on  t hr?  observed  radiation  at 

the  satellite  point.  Thus,  the  upwelliny  radiance  for  the  6-7  pm 

H^O  channel  is  only  affected  by  the  high  cloud.  Consequently,  the 
he 

term  I-  in  £q.  (3.22)  may  be  considered  as  the  satellite  observed 
radiance.  Since  the  cloud  top  height  z^  has  been  estimated  from  the 
foregoing  analysis,  the  cirrus  cloud  thickness  Az^  may  be  subsequently 
determined  using  Eq.  (3.22). 

3.2.3  Determination  of  low  cloud  top  height.  Examining  the  weighting 
function  depicted  in  Fig.  3.2,  we  see  that  most  of  the  energy  contri¬ 
bution  to  the  upwelling  radiance  for  the  11  pin  window  channel  (channel 
8)  comes  from  the  atmosphere  below  700  mb.  If  a  low  cloud  has  an 
optical  depth  large  enough  to  be  considered  a  blackbody  radiator, 
then  the  total  upwelling  radiance  at  the  window  channel  will  arise 
from  two  energy  sources;  one  is  the  gaseous  emission  of  the  clear 
atmosphere  above  the  cloud  top,  while  the  other  is  the  blackbody  emission 
of  the  cloud  layer.  If  the  temperature  profile  and  the  average  cloud 
temperature  are  known,  then  there  should  exist  a  simple  relationship 
between  the  upwelling  radiance  and  cloud  top  height. 

However,  before  finding  the  "observed  radiance",  which  is 
only  affected  by  the  low  cloud,  we  need  to  remove  the  cirrus  cloud 
effect  from  the  actual  observed  radiance.  After  the  cirrus  cloud  top 
height,  cloud  thickness,  and  cloud  properties  at  the  11  inn  window 

R 

channel  have  been  obtained,  the  reconstructed  upwelling  radiance  I.(m) 
may  be  derived  from  Eq.  (3.8)  as  follows: 
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I(  ■  '  '?H  -  >2  -  ['-vf(Azh)].'t('-h)8itT(zh)]',/7f(Azh)  .  I2 


l2  I  Bj[T(z)]Kj(z)dz  ,  (3.?3) 

zh 

where  I ^ (•*> )  is  the  satellite  observed  radiance.  For  the  window 
channel,  we  may  assume  that  the  atmosphere  above  the  cirrus  cloud  top 
is  transparent,  i.e.,  '/’•  (",zh):'l ,  or  K^(z)=0.  So  of  Eq.  (3.23)  is 
approximately  equal  to  zero,  and  Eq.  (3.23)  reduces  to  the  form 


I-  =  I-R  -  [l-^C(Azh)]!Ti (■'‘•zh)]>/y.5C(Azh)  -  (3.24) 

R 

After  reconstructing  the  "observed  radiance",  we  may  then  relate  I- 
and  the  low  cloud  top  height. 

To  a  good  approximation,  water  clouds  may  be  considered  as 
blackbody  radiators.  Thus,  the  cloud  transmissivity  term  in  Eq.  (3.15) 
may  be  ignored  and  we  get 

rr> 

I-  =  B.(Tc)r.(«*,ze)  +  /  Bi[T(z)]K.(z)dz  .  (3.25) 

z,l 

Utilizing  a  prescribed  temperature  profile  and  assuming  that  the 

low  cloud  temperature  Tc  is  the  same  as  the  average  ambient  temperature, 
R 

is  clearly  a  function  of  the  low  cloud  top  height  z? .  On  the  basis 

~  R 

of  the  theoretical  analysis,  the  relationship  between  1^  and  z^  is 
approximately  given  by  the  following  polynomial  equation: 

^  =  e0  +  e,  +  e2Uj)2  ,  (3.26) 
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where  the  coefficients  e^,  ,  and  are  determined  through  empirical 

regression  analyses.  Using  Eq.  (3.26)  the  low  cloud  top  height  may 
be  estimated  from  the  observed  radiance  in  the  window  channel. 

3.3  Parameterization  of  Microwave  Radiative  Transfer 

in  Cloudy  Atmospheres 

In  an  overcast  precipitating  atmosphere  the  cloud  top  height 
may  be  considered  to  be  the  cloud  thickness.  Moreover,  the  cloud  re¬ 
flectivity  may  be  ignored.  Thus,  the  upwelling  brightness  temperature 
in  cloudy  atmospheres  at  the  satellite  point  may  be  decomposed  into 
five  sources  as  schematically  illustrated  in  Fig.  3.3.  These  are  (a) 
the  direct  contribution  from  the  surface,  (b)  emission  from  the  cloud 
layer,  (c)  reflection  of  the  cloud  emission  from  the  surface,  (d)  the 
gaseous  emission  from  the  atmosphere  above  cloud,  and  (e)  reflection  of 
the  gaseous  emission  above  cloud  from  the  surface.  Thus,  the  brightness 
temperature  in  cloudy  atmospheres  may  he  expressed  in  terms  of  the 
parameterized  equation  in  the  form 


W  =  tl>  ..Te  +  +  ^TJr.KzJ 


v>  s 


v  C 


V'  J  V 


VC  V 


OO  (O  lj‘  ?  2  ) 

+  /  T(z)K„(z)dz  +  (2*)2(1-.:v)  /  [ry,-zy  ]2  T(z)Kv(z)dz  , 


Z*.  v 


(3-27) 


where  denotes  the  cloud  emissivity  assuming  a  constant  cloud  tempera¬ 
ture,  ev  is  the  surface  emissivity,  the  cloud  transmissivity,  Tc  the 
average  cloud  temperature,  which  is  assumed  to  be  the  same  as  the  average 
ambient  temperature,  7^  the  gaseous  transmittance,  z^  the  low  cloud 
top  height,  and  K  the  weighting  function,  which  is  in  reference  to 


the  top  of  the  atmosphere. 


Upon  utilizing  a  comprehensive  and  rigorous  microwave  radiative 
transfer  program  developed  at  the  University  of  Utah,  the  theoretically 
simulated  data  may  be  generated  and  regression  eguations,  which  express 
>  c,  TC ,  and  their  product  i  S'c  as  functions  of  the  liquid  water  content 

V  V  V  V 

Aw  and  cloud  thickness  A z,  maybe  constructed.  By  means  of  this  theoreti¬ 
cal-empirical  analysis,  the  following  equations  are  obtained: 

cC  =  a^  +  a^Aw  +  a2AwAz  +  a^Az  +  a^A z)^  ,  (3.28a) 

=  bg  +  b,Aw  +  bpAwAz  +  by\z  +  b^(.Az)^  ,  (3.28b) 

c^rC  =  cQ  +  c^Aw  +  C^AwAz  +  C-jAz  +  c^(Az)?  ,  (3.28c) 

Note  that  the  above  equations  are  linear  with  respect  to  Aw.  The 
computed  emissivity  and  transmissivity  are  highly  correlated  with  the 
real  values.  The  multiple  correlation  coefficient  using  Eq.  (3.28a) 
and  (3.28b)  is  equal  to  0.999,  while  for  Eq.  (3.28c)  it  is  0-969.  The 
latter  correlation  coefficient  is  poorer  than  the  former  because  the 
product  > .  j  v  increases  and  then  decreases  with  increasing  .'w  and  z. 

Since  the  only  term  containing  ;  ^rc  in  Eq.  (3.27)  is  the  surface 
reflection  of  the  cloud  emission  which  contributes  no  more  than  3% 
to  the  total  upwelling  brightness  temperature,  the  error  produced 
in  the  analysis  using  Eq.  (3.28c)  is  insignificant. 

Upon  substituting  Eq.  (3.28)  into  Eq.  (3.27),  calculations 
may  be  carried  out  for  the  brightness  temperature.  Comparisons  with 


values  obtained  from  our  microwave  radiative  transfer  program 
reveal  that  the  accuracy  of  the  brightness  temperature  computed  from 
the  parameterized  equation  is  within  about  ]%.  Moreover,  the  contri¬ 
bution  of  the  last  term  in  Eq.  (3.27)  for  the  water  vapor  and  window 
channels  is  less  than  about  0.1",  and  therefore  it  may  be  neglected  in 
the  numerical  analysis. 

Using  Egs.  ( 3. 28),  emissivity  and  transmissivity  are  plotted 
as  functions  of  the  rainfall  rate  and  cloud  thickness  for  the  SCAMS 
water  vapor  channel  (Fig.  3.4).  The  similar  behavior  for  these  values 

can  be  found  in  the  window  and  oxygen  channels.  In  Fig.  3.4,  it  is 

c  c 

seen  that  7  decreases  with  increasing  Aw  and  Az,  whereas  e  increases 

v  v 

with  increasing  \w  and  Az.  Their  product  (see  Fig.  3.5)  also  increases 

with  increasing  Aw  and  Az,  implying  that  is  a  dominant  factor  when 

the  product  of  the  rainfall  rate  (RFR)  and  thickness  Az  is  less  than 

30.  When  the  product  of  RFR  and  Az  is  larger  than  30,  begins  to 

decrease  with  increasing  Aw  and  Az.  This  implies  that  T ^  is  now  a 

dominant  factor.  If  the  product  of  RFR  and  Az,  that  is  the  cloud 

c  c 

opacity,  becomes  infinite,  clearly  then  the  product  of  T ^  and  r 
equals  zero.  It  should  be  noted  that  rainfall  rate  in  mm/hr  can  be 

3 

converted  to  liquid  water  content  in  g/m  through  the  Marshall -Palmer 
size  distribution.  The  relation  between  rainfall  rate  and  liquid  water 
content  is  tabulated  in  Table  3.1. 

3 . 4  Determination  of  L iqu id  Water  Content  from  SCAMS  Channel s 

The  cloud  liquid  water  content  information  is  contained  in 
Eq.  (3.27)  through  values  of  iC  and  7'c,  which  in  turn  are  functions  of 

V  V 
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CLOUD  THICKNESS  (km! 


CLOUD  THICKNESS  (km) 


Fig.  3.4. 


The  cloud  transmissivity  arid  emissivity  as  Functions  of  the 
rainfall  rate  and  cloud  thic knevc  for  the  FOAMS  water  vapor 
c ha nne  I . 


CLOUD  THICKNESS  (km) 


The  product  of  the  cloud  transmissivity  and  emissivity  as  a 
function  of  the  rainfall  rate  and  cloud  thickness  for  the 
CAMS  -water  vapor  channel. 


Table  3.1.  Relation  between  rainfall  rate  and  liquid 
water  content. 


Rainfall  Rate 
(mm/hr) 

1 

2 

3 

r 

4 

5 

6 

7 

8 

9 

10 

Aw(xlO-1 ) 
(g/m3) 

h  H 

0.889 

1.592 

L . 

2.238 

2.849 

3.437 

4 . 006 

4.559 

5.101 

_ 

5.631 

6.152 

i 

_ l 

Aw  and  Az  as  indicated  in  Eq.  (3.28).  Substituting  Eq.  (3.28)  into 
Eq.  (3.27)  and  after  a  number  of  mathematical  manipulations,  Eq.  (3.27) 
becomes 


Tg(v)  =  [u]  +  u2(zy)]Aw  +  uQ  +  u3(z?)  + 


u4^zi ^ 


(3.29) 


where 


Tg(v)  =  [Tg(v)  ‘  /  T(z)Kv(z)dz]/7’vKzy)  , 

z?, 

and 


ui 


Vc 


biV 


+  ciVWv>’  i=0>1’2*3’4 


If  the  temperature  and  humidity  profiles  are  prescribed,  and  the  cloud 
thickness  or  cloud  top  height  z„  has  been  determined  from  previous 
analyses,  then  Tg(v)  may  be  estimated  from  the  observed  brightness 
temperature  Tg(v).  Furthermore,  if  we  estimate  e  from  the  SCAMS  data 
for  the  water  vapor  and  window  channels  utilizing  an  empirical  equation 
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proposed  by  Waters  et  al.  (1975),  the  liquid  water  content  can  then  be 
obtained  from  Eq.  (3.29)  in  the  form 


Aw  =  Tg( v)  -  l[uQ  +  u3(z{)  +  u4(zv)2]!/[u1  +  u2(z^)]  .  (3.30) 

We  shall  now  describe  the  procedure  for  determining  the  sur¬ 
face  emissivity  from  the  SCAMS  channels.  Since  the  instrumental  char¬ 
acteristics  for  the  Nimbus  6  SCAMS  channels  1  and  2  (at  22.235  and 
31.650  GHz)  are  similar  to  those  for  the  Nimbus  5  NEMS  channels  1  and 
2  (at  22.235  and  31.4  GHz),  the  equation  developed  by  Waters  et  al. 
(1975)  for  calculating  surface  reflectivities  may  be  used  in  the 
present  numerical  analysis.  The  surface  emissivity  over  land  can  be 
simply  estimated  from  the  surface  reflectivity  R  in  the  form 

=  1  -  Ry  .  (3.31) 


According  to  Waters  et  al.,  the  surface  ref  1  ect’vi  ty  forchannel  2  can 
be  approximated  by 


[Tb(2)  -  Te(2)]  +  [Te(2)  -  TJ]  exp  (-i°) 
' Tg( 2) [ 1 -exp( - i ^ ) 1  '  V  exp(-.°) 


(3.32) 


where  Tp(2)  is  the  measured  brightness  temperature  for  channel  2.  Tg(2) 

is  an  equivalent  atmospheric  emission  temperature  for  channel  2,  which 
is  approximately  given  by 


T  (2)  ;  0.965  x  T(850  mb) , 


(3.33) 
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where  T { 850  mb)  is  the  prescribed  atmospheric  temperature  at  the  850  mb 
pressure  level.  The  atmospheric  transmissivity  for  channel  2,  exp  (-ifj), 
is  determined  from  the  empirical  equation 


ex  p ( - ' p) 


b  ' l-a[r 


T.s-V1)] 

Ts  -  hW 


M2 


(3.34) 


where  a=0.32,  b=0.66,  and  ,  R^/Rj ,  which  has  a  value  of  0.94  and  1.0 
for  flat  water  surface  and  land,  respectively,  and  Tg(l)  is  the  mea¬ 
sured  brightness  temperature  for  channel  1.  By  virtue  of  Eqs.  (3.32) 
-  (3.34),  the  surface  reflectivity  for  channel  2  can  be  computed 
and  the  surface  emissivity  thus  obtained  from  Eq.  (3.31). 


3 ^  Resu 1 ts  of  Theoretical  Ana  1 yses 

In  order  to  test  the  validity  of  the  theoretical  analyses 
outlined  previously,  a  hypothetical  atmosphere  containing  two-layered 
clouds  is  constructed.  The  cirrus  cloud  top  is  assumed  to  be  located 
at  300  mb  (9.358  km)  for  the  mid-latitude  summer  profile  with  a  thickness 

of  1  km.  The  low  cloud  top  height  is  placed  at  780  mb  haviny  a  total  liquid 

2 

water  content  of  0.07  g/cm  .  The  "observed"  infrared  upwelling  radiances  for  HIRS 
channels  are  computed  from  Eq.  (3.14),  while  the  "observed"  microwave 
upwelling  brightness  temperatures  for  SCAMS  channels  are  calculated 
from  Eq.  (3.27).  Upon  imposing  different  percentages  of  random  errors 
on  the  "observed"  values,  computations  of  the  cloud  parameters  are 
subsequently  carried  out.  Results  of  the  sequential  calculations  for 
the  cirrus  cloud  top  height,  cloud  thickness,  low  cloud  top  height,  and 
cloud  liquid  water  content  are  shown  in  Figs.  3.6.  For  the  random 
error  ranging  from  0  to  37.,  the  range  of  the  cirrus  top  height  is 
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from  9-6  to  7.0  km.  The  computed  height  errors  increase  linearly  with 
increasing  random  error  as  shown  in  Fig.  3.6a.  Note  that  without  impos¬ 
ing  random  errors  on  the  "observed"  values,  the  computed  cirrus  top 
height  with  a  value  of  9.6  km  is  slightly  higher  than  the  assumed  value. 
This  is  due  to  interference  from  low  clouds  on  channels  4  and  5  of  the 
HIRS,  which  are  used  to  estimate  the  cirrus  top  height.  But  this  slight 
difference  is  not  significant  enough  to  cause  errors  on  the  sequential 
calculations  in  the  numerical  experiment.  The  increase  of  the  competed 
errors  for  the  cirrus  thickness  is  also  linear  with  respect  to  the  in¬ 
crease  of  the  imposed  random  error.  However,  the  error  range  of  the 
thickness  is  limited  to  within  about  0-07  km  with  an  assumed  thickness 
of  1  km.  The  low  cloud  top  height  is  calculated  after  the  cirrus  cloud 
effect  has  been  removed  from  the  HIRS  window  channel  measurement.  As 
shown  in  Fig.  3.6b,  the  low  cloud  top  height  increases  when  the  cirrus 
top  height  and  thickness  decrease  due  to  the  increased  random  error. 

This  may  result  from  compensating  effects  in  which  the  cirrus  parameters 
are  underestimated  while  an  overestimation  of  the  low  cloud  thickness 
produces  the  same  upwelling  radiance.  The  dashed  line  in  the  figure 
shows  the  computed  liquid  water  content.  The  fluctuation  of  the  liquid 
water  content  reveals  a  nonlinear  dependence  of  the  liquid  water  con¬ 
tent  on  other  parameters.  As  can  be  seen,  even  without  inserting  ran- 
dom  errors,  the  computed  liquid  water  content  (0.073  g/cm  differs  from 
the  assumed  value  of  0.07  g/cm2.  However,  when  random  errors  are  inserted 
‘  '>■  computed  liquid  water  content  remains  in  the  range  of  0.08  to  0.066 
This  range  appears  to  be  acceptable  in  view  of  the  difficulty 
.  deriving  the  liquid  water  content  from  any  remote 
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sensing  and  in-situ  technique. 


3 . 6  Appl i cat ions  to  Nimbus  VI  I1JRS and  SCAMS  Data 

In  this  section,  we  describe  the  resulting  cloud  parameters 
inferred  from  Nimbus  VI  HIRS  and  SCAMS  data  for  two  different  weather 
systems  documented  in  a  previous  scientific  report  by  Astling  and  Liou 
(1979).  These  weather  systems  were  associated  with  summertime  convec¬ 
tion  and  wintertime  mid-latitude  cyclones  which  occurred  in  the  United 
States  during  August  25,  1979  and  February  20,  1976,  respectively.  On 
August  25,  only  seven  cases  met  the  criteria  that  infrared  and  micro- 
wave  data  from  the  Nimbus  VI  satellite  are  available  along  with  conven¬ 
tional  meteorological  observations.  On  February  20,  1976,  six  cases 
were  selected. 

The  cloud  thickness  and  liquid  water  content  derived  from 
HIRS  and  SCAMS  data  for  the  August  25  cases  are  listed  in  Table  3.2. 

On  this  date,  the  major  surface  weather  features,  the  500  mb  flow 
patterns,  the  GOES  visible  and  infrared  imagery,  the  radar  sunmaries, 
and  the  precipitation  amounts  for  six  hour  periods  were  comprehensively 
described  in  the  report  by  Astling  and  Liou. 

Surface  reports  indicated  the  rainfall  rate  at  Lake  Charles, 

Louisiana,  was  about  1  mm/hr  at  1800  Z.  The  computed  liquid  water 

content,  high  cloud  thickness,  and  low  cloud  thickness  using  satellite 

-2 

data  at  1717  Z  is  0.135  g  cm  ,  0.006  km  and  3.77  km,  respectively. 

The  computed  liquid  water  content,  when  converted  to  the  rainfall  rate 
utilizing  the  Marsha  1 1 -Pa liner  size  distribution,  is  equivalent  to 
5  mm/hr,  which  is  much  larger  than  the  rainfall  rate  observed  at  the 
surface.  Since  .nr  .<e,  radar,  and  synoptic  reports  indicate  strong 
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Table  3.2.  Cloud  thickness  and  liquid  water  content 

derived  from  Nimbus  6  HIRS  and  SCAMS  channel 


STATION 

! 

TIME 

(GMT) 

i 

SURFACE 

REPORT 

i 

i 

CLOUD 

THICKNESS  (km) 
low  high 

i.  .  .  .  _ 

LWC  2 
(g  crn  ) 

Lake  Charles, 
LA 

17:17 
Aug.  25 

Precip.  (1  mm/ hr)  ^ 
Cb  (9/10) 

3.77 

0.006 

0.135 

Peoria,  IL 

17:19 

I  Aug.  25 

Precip.  (0.5  mm/ hr) 
Low  Cloud  (10/10) 

4.07 

0.271 

1 

0.060 

Internat1 1 
Falls,  MN 

17:22 
Aug.  25 

Precip.  (trace)  1 
Low  Cloud/Cirrus 

_ 41 0/10). _ j 

6.48 

0-194 

0.005 

Topeka,  KS 

17:19 
Aug.  25 

1 

Low  Cloud/Cirrus 
(10/10) 

3-42 

0-102 

0.013 

Green  Bay, 

MI 

.  .  1 

17:20 
Aug.  25 

Total  Overcast 
(10/10) 

4.13 

0.336 

0.056 

Dayton,  OH 

17:18 
Aug.  25 

L- 

Low  Cloud 
(10/10) 

3.437 

1 

! 

u 

0.043 

0.037 

Huntington, 

WV 

! 

17:18 
!  Aug.  25 

i 

;  Low  Cloud 

1  (10/10) 

i 

3.18 

0.055 

0.016 

convective  activity  in  the  area  the  observed  cumulonimbus  clouds  could 
have  contained  larger  cloud  droplets.  Thus,  it  appears  that  the  larger 
liquid  water  content  derived  from  the  microwave  channels  may  well 

represent  the  actual  moisture  present  within  the  cumulonimbus. 

_2 

We  found  a  liquid  water  content  value  of  0.005  g  cm  for 
the  International  falls  case.  Our  computed  value  matches  quite  well 
with  the  surface  precipitation  report.  Although  the  low  cloud  thick¬ 
ness  is  estimated  to  be  about  6.5  km,  the  cloud  might  well  contain  a 
large  number  of  ice  crystals  and  small  water  droplets.  Under  this 
condition  only  large  raindrops  near  the  cloud  bottom  would  contribute 
significantly  to  the  signature  in  the  microwave  region. 

At  1719  Z  25  August,  according  to  an  1800  Z  surface  report 
and  1900  Z  GOES  infrared  imagery,  Topeka,  Kansas  was  covered  by  a  low 
cloud  deck  with  thin  cirrus  aloft.  Inspection  of  the  sequential 
surface  maps  reveals  a  major  cold  front  passed  Topeka  at  about  1200  Z. 
As  described  in  the  report  by  Astling  and  Liou,  clouds  in  the  post 
frontal  region  tend  to  be  less  convective  and  contain  smaller  amounts 
of  moisture.  Thus,  the  computed  cloud  liquid  water  content  of  0.013  g 
cm  (equivalent  to  0.4  mm/hr  rainfall  rate)  with  3.42  and  0.1  km  for 
low  and  high  cloud  thickness,  respectively,  appears  to  agree  well 
with  synoptic  observations. 

The  precipi tation  report  from  Peoria,  Illinois  at  1800  Z  was 

0.5  mm/hr  rainfall  rate.  The  calculated  liquid  water  content  is 
-2 

0.06  g  cm  ,  which  is  equivalent  to  less  than  1  mm/hr  rainfall  rate. 
Again,  the  larger  computed  value  for  the  cloud  liquid  water  content 
may  be  associated  with  a  wel 1 -developed  convective  cloud. 


The  strong  frontal  activity  did  not  produce  precipitation  at 

Green  Bay,  Wisconsin  or  Dayton,  Ohio.  However,  it  is  believed  that  a 

large  amount  of  moisture  was  present  in  the  atmosphere  in  the  Great 

Lakes  area.  The  computed  cloud  liquid  water  content  for  the  Green 

_2 

Bay  case  is  0.056  g  cm  ,  which  would  produce  1.7  mm/hr  rainfall  rate 

if  the  total  cloud  moisture  were  converted  to  precipitation.  The 

_2 

cloud  liquid  water  content  for  the  Dayton  case  is  0.037  g  cm  , 
equivalent  to  1.3  mm/hr  rainfall  rate. 

The  overcast  low  cloud  reported  at  Huntington,  West  Virginia, 
did  not  contain  as  much  moisture  as  the  Great  Lakes  cases.  The  cal¬ 
culated  values  show  the  cloud  liquid  water  content  at  Huntington 

_2 

was  about  0.016  g  cm  ,  which  is  approximately  equal  to  a  0.5  mm/hr 
rainfall  rate.  The  small  cloud  liquid  water  content  obtained  for  this 
case  might  be  due  to  the  lack  of  a  moisture  supply,  since  the  station 
was  far  east  of  the  very  active  frontal  zone. 

The  cloud  thickness  and  liquid  water  content  derived  from 
HIRS  and  SCAMS  data  for  the  February  20  and  21  cases  are  listed  in 
Table  3.3.  In  reference  to  the  synoptic  analyses  discribed  by  Astling 
and  Liou,  the  major  trough  with  low  center  located  in  Kansas  brought 
abundant  amount  of  moisture  from  the  Gulf  of  Mexico  to  Lousiana  and 
into  Arkansas  through  the  warm  front.  The  surface  report  of  Little 
Rock,  Arkansas  at  1800  Z  20  February,  indicated  a  total  obstruction 
due  to  low  cloud.  The  computed  results  revealed  that  the  cloud  de¬ 
velopment  at  Little  Rock  was  quite  stable  since  the  low  cloud  thickness 

only  extended  to  about  1.46  km.  The  total  cloud  liquid  water  content 

_2 

was  estimated  to  be  about  0.0067  g  cm  ,  which  is  equivalent  to  about 
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Table  3.3 


Cloud  thickness  and 
derived  from  Nimbus 
for  February  20  and 


1 iquid  water  content 
VI  HIRS  and  SCAMS  channels 
21  cases. 


STATION 

TIME 

SURFACE  REPORT 

CLOUD 

LWC 

(GMT) 

r 

I 

THICKNESS  (km) 

(g  cm“2) 

J 

1 

i 

Low  cloud  (10/10) 

low  high 

Little  Rockl 
ARK  ; 

17:11 
Feb.  20 

1.46  0.05 

0-0067 

Rapid  City, 
S.D. 

17:14 
Feb.  20 

i  Low  cloud  (10/10) 

1  1 

3.23  0.47  ' 

0.01 

- . — - 

- - - - - | 

— - 

— 

Spokane 

10:02 

Low  cloud  (10/10) 

5.09  0.03 

0.22 

WA. 

Feb.  20 

with  fog 

Huntington, 
W.  VA. 

16:30 
Feb.  21 

Low  cloud  (10/10) 

1.17  0.92 

0.078 

Dayton, 

OH 

16:30 
Feb.  21 

Precipitation  (2  mm/hr) 
Low  cloud  (10/10) 

3.85  0.39 

0.30 

Peoria,  IL. 

16:31 
Feb.  21 

Precipitation  (0.5  mm/hr) 
Low  cloud 

3.88  0.52 

0.40 
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0.7  mm/hr  rainfall  rate. 

At  the  same  time,  a  moist  tongue  at  850  mb  stretched  from 

Gulf  coast  along  eastern  Texas  and  Lousiana,  across  Arkansas,  Iowa, 

into  southern  South  Dakota.  Rapid  CHy,  South  Dakota,  reported  an 

overcast  low  fractostratus  and/or  fractocumulus  condition  with  a 

ceiling  near  1,000  ft.  according  to  the  1800  Z  surface  report.  These 

cloud  types  are  associated  with  the  approaching  cold  front.  The 

results  we  obtained  showed  that  Rapid  City  was  covered  by  a  low  cloud 

deck  having  a  3.2  km  thickness  and  a  thin  cirrus  cloud  with  a  0.47  km 

thickness.  The  cloud  liquid  water  content  was  found  to  be  about 
_2 

0.01  g  cm  ,  which  is  equivalent  to  about  0.35  mm/hr  rainfall  rate. 

In  addition  to  the  low  pressure  centered  in  Kansas,  another 
low  pressure  was  located  just  off  the  northeast  Pacific  Ocean.  This 
system  moved  southeastward  into  the  state  of  Washington.  The  strong 
pressure  gradients  and  cyclonically  curved  isobars  brought  a  signifi¬ 
cant  amount  of  Pacific  moisture  ashore.  The  surface  weather  observa¬ 
tion  at  Spokane,  Washington  at  1800  Z  20  February  reported  a  low  cloud 
total  obstruction  with  cloud  ceiling  height  at  about  300  ft.  Ground 
fog  and  near  zero  visibility  were  also  reported.  The  computed  results 

indicated  that  a  low  cloud  with  a  thickness  of  5.09  km  and  liquid 

_2 

water  content  of  0.22  g  cm  was  present  at  Spokane  under  the  satellite 
pass  at  1902  Z. 

The  precipitation  activity  was  most  intense  along  the  occluded 
and  cold  fronts  located  in  the  midwest  during  the  12  hours  prior  to 
0000  Z  22  February.  Six  hour  rainfall  amounts  were  found  to  be  in 
excess  of  1  inch  at  several  stations  along  the  cold  front. 


At  Dayton,  Ohio,  a  2  inm/hr  rainfall  rate  and  a  low  cloud 
ceiling  near  400  ft.  were  reported  at  1800  Z  21  February.  The  com¬ 
puted  results  show  that  there  were  thick  low  clouds  having  total  liquid 
water  content  of  as  much  as  0.3  g  cm  (equivalent  to  about  10  mm/hr 
rainfall  rate).  The  results  appear  to  be  in  qualitative  agreement  with 
the  surface  and  synoptic  reports  which  indicated  strong  ascending 
motions  at  Dayton. 

The  retrieved  results  for  Peoria,  Illinois  show  a  low  cloud 

_2 

deck  3.88  km  thick  with  a  large  liquid  water  content  of  0.4  g  cm 
(11  mm/hr  rainfall  rate).  This  value  is  much  larger  than  the  observed 
rainfall  rate  of  0.5  mm/hr.  However,  according  to  the  index  of  stability 
from  the  synoptic  data,  the  area  showed  strong  ascending  motions  which 
could  produce  large  water  droplets  within  the  cloud,  and  hence,  large 
liquid  water  content. 

The  retrieved  cloud  parameters  for  the  Huntington,  Virginia 

case  reveal  that  the  low  and  high  cloud  thicknesses  are  1.17  and  0-92  km, 
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respectively,  and  the  liquid  water  content  is  0.078  g  cm  .  The  reason 
for  the  small  liquid  water  content  is  probably  due  to  the  relatively 
stable  atmosphere  surrounding  the  area.  The  radiosonde  observations 
indicated  that  the  relative  humidity  was  90%  at  300  mb  but  was  less 
than  50%  at  850  mb. 

In  summary,  we  have  applied  the  cloud  sounding  technique 
developed  from  the  parameterization  of  infrared  and  microwave  radiative 
transfer  to  the  Nimbus  VI  DIRS  and  SCAMS  data.  The  satellite  retrieved 
cloud  parameters  from  selected  infrared  and  microwave  channels  appear 
to  be  qualitative  agreement  with  the  available  synoptic,  radiosonde. 


surface  and  radar  reports.  The  results  illustrated  in  this  section 
are  indeed  encouraging.  Thus,  it  appears  that  applications  of  the 
cloud  sounding  technique  developed  in  this  report  for  the  mapping  of 
three-dimensional  cloud  compositions  and  structure  over  the  globe, 
utilizing  the  available  satellite  infrared  and  microwave  data,  may 


be  carried  out. 


SECTION  4 


CONCLUSIONS 

In  this  report,  we  have  presented  two  areas  of  the  research 
results  concerning  the  sounding  of  cloud  parameters  from  orbiting 
meteorological  satellites.  First,  in  conjunction  with  the  cloud  sounding 
utilizing  the  Nimbus  6  ESMR  channel,  we  have  developed  a  microwave  radia¬ 
tive  transfer  program  in  which  the  polarization  effect  is  included  so 
that  the  vertically  and  horizontally  polarized  brightness  temperatures 
can  be  computed  and  compared  with  observed  data.  The  program  involves 
the  use  of  the  di screte-ordi nate  method  in  solving  the  fundamental  trans¬ 
fer  equation  governing  absorbing  molecules  and  Rayleigh  particulates  in 
the  microwave  region.  Analytical  solutions,  assuming  an  isothermal  cloud 
temperature,  are  derived  for  both  vertical  and  horizontal  polarization 
components  in  which  the  polarization  cnaracteri sti cs  of  the  surface 
emissivity  is  taken  into  accori.. 

Con  par i sons  between  the  theoretically  computed  brightness 
temperature,  employing  the  available  tanperature  and  mixing  ratio  profiles 
from  radiosonde,  and  the  observed  data  are  carried  out  for  selected  clear, 
cloudy,  and  precipitating  atmospheres  over  inpolarized  land  and  polarized 
ocean  surfaces.  Fairly  good  agreement  is  found  between  the  computed 
results  and  observed  data  with  differences  within  about  6°K.  It  may  be 
concluded  therefore  that  the  radiative  transfer  program  developed  in  this 
study  appears  to  be  quite  reliable  despite  the  difficulty  involved  in  the 
determination  of  the  surface  emissivity  in  all  weather  conditions.  The 
microwave  radiative  transfer  program  is  then  used  to  perform  theoretical 
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sensitivity  analyses  to  investigate  the  effects  of  the  emergent  angle, 
rainfall  rate,  cloud  thickness,  and  atmospheric  profile  on  the  brightness 
temperatures  for  vertical  and  horizontal  polarization  components.  Sub¬ 
sequently,  we  develop  a  simple  empirical -theoretical  approach  to  recover 
the  cloud  liquid  water  content  and  cloud  thickness  from  vertical  and 
horizontal  brightness  temperatures.  Using  the  available  ESMR  data  over 
ocean  areas,  we  illustrate  that  reasonable  cloud  liquid  water  contents 
and  cloud  thicknesses  may  be  derived  from  a  set  of  linear  equations. 

However,  owing  to  the  unavailability  of  surface  weather  observations  over 
ocean,  we  have  been  unable  to  verify  the  degree  of  reliability  for  these 
retrieval  analyses. 

It  should  be  pointed  out  here  that  the  ESMR  program  utilizes 
highly  simplified  linear  equations  in  the  retrieval  of  cloud  parameters. 

The  next  hierarchy  of  the  theoreti cal -emf i ri cal  approach  would  be  the 
development  of  high  order  equations  for  retrieval  purposes  on  the  basis 
of  the  physical  and  mathematical  analyses  of  the  solution  of  the  microwave 
radiative  transfer  equation.  Moreover,  one  additional  improvement  to  the 
theoretical  microwave  radiative  transfer  program  would  be  to  remove  the 
Rayleigh  approximation  in  the  treatment  of  the  source  function.  It  is 
believed  that  the  analysis  involving  the  Nimbus  6  ESMR  data  illustrated  in 
this  report  may  be  of  direct  use  to  the  DMSP  SSM/I  multi  spectral  microwave 
imager  to  be  flown  in  1983  for  applications  and  operational  uses. 

The  second  area  of  research  results  presented  in  this  report  is 
concerned  with  the  exploitation  of  cloud  sounding  techniques  utilizing  a 
combination  of  infrared  and  microwave  channels.  In  connection  with  this 
purpose,  we  have  first  developed  parameterized  equations  for  infrared  and 
microwave  radiation  transfc  in  atmospheres  containing  multi-layered  clouds. 


87 


Parameter-! zations  of  the  infrared  and  microwave  radiation  transfer  makes 
use  of  the  physical  and  geometrical  principles,  and  they  are  verified 
through  rigorous  transfer  programs  developed  at  the  University  of  Utah 
in  the  last  few  years.  The  parameterized  equations  are  then  employed  to 
retrieve  high  cloud  top  height,  high  cloud  thickness,  low  cloud  top 
height,  and  cloud  liquid  water  content  successively  utilizing  selected 
HIRS  and  SCAMS  channels. 

Hypothetical  error  analyses  have  been  performed  and  we  find 
that  the  retrieval  technique  involving  the  use  of  both  infrared  and  micro- 
wave  channels  appears  to  be  theoretically  rigorous  and  practically  workable. 
Consequently,  we  have  applied  this  cloud  sounding  technique  to  the  Nimbus  6 
HIRS  and  SCAMS  data  for  a  number  of  carefully  selected  cases  associated 
with  summertime  convective  cloud  systems  and  wintertime  large  scale  synoptic 
cyclones  described  in  the  report  by  Astling  and  Liou  (1979).  High  and  low 
cloud  thicknesses  and  total  liquid  water  content  are  recovered  from  a 
combination  of  selected  HIRS  and  SCAMS  channels  and  are  compared  with  the 
available  synoptic,  radiosonde,  surface,  and  radar  observations  and  satel¬ 
lite  photographs.  Qualitative  agreement  is  found  between  the  retrieved 
cloud  parameters  and  the  results  estimated  from  synoptic  data  and  analyses. 

The  development  of  the  cloud  sounding  technique  using  combined 
infrared  and  microwave  channels  is  quite  significant  in  view  of  the  rela¬ 
tive  advantages  and  shortcomings  of  the  infrared  and  microwave  frequencies 
with  respect  to  cloud  sounding  problem.  It  may  be  pertinent  to  point  out 
here  that  microwave  frequencies  penetrate  through  low  clouds,  which  are 
opaque  in  the  infrared  region,  and  that  high  clouds  normally  affect 
significantly  the  transfer  of  infrared  radiation  but  have  little  effect  on 
the  microwave  propagation.  It  is  quite  clear  therefore  that  a  combination 
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of  infrared  and  microwave  frequencies  is  vitally  essential  to  the  determi¬ 
nation  of  the  vertical  cloud  structure  and  compositions.  In  this  investiga¬ 
tion,  we  have  demonstrated  convincingly  the  possibility  of  deriving  the 
vertical  cloud  parameters  based  on  a  number  of  carefully  selected  cases. 
Thus,  further  investigations  and  verifications  may  pave  the  way  for  the 
mapping  of  three-dimensional  cloud  parameters  over  the  global  space. 

Finally,  it  should  be  pointed  out  that  the  fundamental  difficulty  in  the 
development  of  the  cloud  sounding  technique  from  orbiting  meteorological 
satellites  has  been  the  lack  of  cloud  physics  data  for  intercomparisons 
and  verifications.  In  this  connection,  it  appears  logical  to  suggest 
that  a  carefully  designed  field  experiment  involving  the  simultaneous 
satellite  radiation  program,  aircraft  cloud  physics  program,  and  synoptic 
and  radiosonde  observations  be  carried  out  when  the  developed  cloud  sound¬ 
ing  technique  using  both  infrared  and  microwave  sounders  appears  promising. 
The  accurate  and  reliable  cloud  sounding  program  in  an  operational  mode 
certainly  will  enhance  and  improve  the  ongoing  AFGWC  3D  NEPH  program. 
Moreover,  it  also  is  vitally  important  to  the  understanding  of  the  radia¬ 
tive  budget  and  climate  of  the  earth-atmosphere  system. 
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